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The earliest attempt to subject the Theory of the Tides to a rigorous dynamical 
treatment was given by Laplace in the first and fourth books of the 'Mecanique 
Celeste/ The subject has since been treated by Airy, # Keltin/I' Darwin,} Lamb,§ 
and other writers, but with the exception of the extension of Laplace's results to 
include the theory of the long-period tides, but little practical advance has been made 
with the subject, in spite of the enormous increase in the power of the mathematical 
resources at our disposal, and the problem has remained in very much the same 
condition as it was left by Laplace. This arises no doubt partly from the difficulties 
inherent to the subject, but partly from the form in which the theory was originally 
presented by Laplace in the f Mecanique Celeste/ which has been described by Airy 
as "perhaps on the whole more obscure than any other part of the same extent in 
that work/' The obscurity complained of does not however seem to have been 
entirely removed by Laplace's successors, and it was the fact that every presentment 
of the theory with which I was acquainted offered some points of difficulty, that in 
the first instance led me to take up the problem ab initio, partly with the purpose of 
allaying the doubts which had arisen in my own mind as to the validity of certain 
approximations employed by Laplace and adopted by his successors, and partly in 
the hope that I might be able to extend the results of Laplace to meet more fully 
the case presented by the circumstances actually existent in Nature. 

Up to the present I have been unable to free the problem to any extent from the 
limitations which have been imposed by previous writers, and consequently it would 
be futile to claim that the results I am now able to put forward materially advance 

* ( Eneyc. Mefcropolitana ' ; Art., " Tides and Waves," Section III. 

t c PM. Mag.,' 1875, vol. 50. 

% * Encyc. Britannica ' (9th edition) ; Art., u Tides." 

§ ' Hydrodynamics,' chapter viii. 
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our knowledge of the tides as they actually exist; but I venture to hope that these 
results, as applied to the oscillations of an ideal ocean, considerably simpler in 
character than the actual ocean, may prove of some interest from the point of view 
of pure hydrodynamical theory. 

In §§ 1-4 I have devoted considerable space to the formation of the dynamical 
equations. The equations obtained agree with those used by Laplaoe, and 
consequently it may be thought that I have been unnecessarily diffuse over this 
part of the subject. My apology is that the questionable, if not erroneous, reasoning 
which has often been assigned for the various approximators introduced seemed to 
me to warrant a very minute examination of the formation of these equations. An 
analytical treatment, such as that I have used, seems to me to be the only safe 
method of procedure to ensure that the approximations do not involve the neglect 
of any terms which may be of equal importance with those retained, many of which 
are extremely small. The method adopted follows Laplace in so far as it consists 
of a transformation of the general equations of oscillation of a rotating fluid, I 
trust however that these general equations in the form I have used, which seems 
to be the simplest form to which they can be reduced, may be found less "repulsive" 
than those employed by Laplace. 

"§ 5 deals with the integration of these equations. The forms of solution discussed 
in the present paper are those which are symmetrical with respect to the axis of 
rotation. The types of oscillation represented by these solutions have been named 
by Laplace, " Oscillations of the First Species," but he omitted to discuss them in 
detail on the grounds that the oscillations of such character, which might be 
expected to exist in Nature, would be modified to such an extent by friction that 
they would be far better represented by the old " equilibrium theory," than by a 
dynamical theory which failed to take due account of the action of friction. The 
tides in question will be of long period, the shortest of the periods being half a lunar 
month in duration, but Professor Darwin was, I believe, the first to call attention 
to the fact that this length of period will hardly be sufficiently great to render the 
effects of friction of such paramount importance, and hence he added to the work of 
Laplace a discussion of the long-period tides when not subject to frictional 
influences. I have recently attempted to estimate the effects of friction on the tidal 
oscillations of the ocean, # and the results at which I have arrived fully confirm 
the view of Professor Darwin as to the small influence of friction on the lunar- 
fortnightly tides, and render it highly probable that the effects will be almost 
equally slight on the solar long-period tides. 

The method of integration I have employed differs from that used by Darwin, my 
aim having been to express the results by means of series of zonal harmonics instead 
of by the power-series obtained by him. The advantages of this are two-fold; firstly, 

* In a paper read before the London Math. Soc, December 10th, 1896. 
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it allows of our including in our analysis the effects due to the gravitational attrac- 
tion of the water ; and secondly, the convergence of the series obtained will be much 
more rapid. The latter circumstance is of particular value, as it has enabled me to 
treat with considerable success the problem of the free oscillations of the ocean. 

The general problem of the small oscillations of a rotating system possessing a 
finite number of degrees of freedom has been discussed by Thomson and Tait ; # but 
the extension to meet the case where the number of degrees of freedom is infinite 
involves analytical considerations of some delicacy. As a rule, the transition from 
the case of a system with finite freedom to that of a system with infinite freedom is 
effected by the employment of "normal coordinates, "t and the chief difficulty in the 
solution of problems relating to the vibrations of the latter class of system consists in 
the discovery of these coordinates. The researches of Thomson and Tait just men- 
tioned shew however that in a rotating system these normal coordinates do not 
exist, and hence that the methods ordinarily employed to deal with the oscillations 
of a system about a state of equilibrium will no longer suffice for the treatment of 
our problem. In most "gyrostatic" problems which have been solved hithertoj the 
solution has been obtained by means of a system of quasi- normal coordinates. When 
such coordinates exist, only a finite number of oscillations of certain particular types 
are possible, and, by constraining the system to vibrate in one of these types, we may 
treat it in the same manner as a system with a finite number of degrees of freedom. 
The period-equation for the free oscillations of an assumed type will then only possess 
a finite number of roots, and will consequently be an algebraic equation usually 
most readily obtained in a determinantal form. It is shewn at the end of § 5 that 
the coordinates we have used possess this property when the depth of the ocean 
follows certain restricted laws ; but in general no such quasi-normal coordinates exist, 
and whatever coordinates be employed, the displacements in any of the fundamental 
modes of vibration can only be expressed, by means of an infinite number of such 
coordinates. The most advantageous choice of coordinates will then be that which 
leads to most rapidly converging series. 

As however the oscillations of an assumed type can only be expressed by an infinite 
series of coordinates, it follows that an infinite number of oscillations of any assumed 
type must be possible, and that consequently the period-equation for oscillations of 
this type will have an infinite number of roots and will therefore be transcendental 
instead of algebraic in character. It is possible that the transition from systems 
with finite freedom to systems with infinite freedom may be treated with advantage 
by the employment of determinants of infinite order (as a means of expressing the 
transcendental period-equation), after the manner introduced into analysis by 

* ' Natural Philosophy,' Part I., § 345. 

f JUyleigh, ' Theory of Sound,' vol. 1, § 87. 

J Gf. Poincare, ' Acta Mathematical vol. 7 ; Bryan, ' Phil. Trans.,' 1889. 
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G. W, Hill, in his ' Researches on the Lunar Theory ;' # but in the present instance 
we are able to avoid the difficulties involved in the use of these infinite determinants, 
in that the forms of determinant which occur are those which are associated with 
continued fractions. 

§ 6 deals with the analytical discussion of the deduction of the period-equation. 
The method is based on a paper by Lord Kelvin,! in which the author defends the 
procedure of Laplace against certain allegations to which it had been subjected 
by Airy, but I have endeavoured to present the arguments in a somewhat different 
light, so as to bring out more clearly the analogy between our problem and the 
general problem of vibrating systems with finite freedom. 

Tn §§ 7-10 I have given illustrations of the method of solving the period-equation 
numerically, and of the subsequent determination of the type of motion for the 
different fundamental modes. As the ground covered in these sections is almost 
entirely new, I have devoted considerable time and labour to the arithmetical 
determination of the periods and types of the principal oscillations for a system 
comparable with the earth in magnitude. The results are tabulated in these 
sections. 

§11 deals briefly with the forced tides of long period due to the moon in an ocean 
of uniform depth. The results agree with those previously obtained by other 
methods, but differ from them in analytical form. In § 12 I have given illustrations 
of a means of extending the method of numerical computation to cases where the 
law of depth is less restricted in character. 

The consideration of forced tides of very long period, dealt with in § 13, points to 
the existence of free oscillations of infinitely long period. This, I believe, was 
first noted by Professor Lamb, J but the application of the dynamical equations for 
the tides to the treatment of these free oscillations has not been previously carried 
out. The types of motion in question appear to be of considerable importance, as 
they throw light on a phenomenon which in the past has been the subject of 
considerable controversy. The difficulties which have been met with in attempts to 
account for the existence of ocean currents all seem to me to arise from an over- 
estimate of the effects of viscosity on the motion of the sea. The large-scale ocean 
currents have been attributed by Sir John Herschel§ and others entirely to the 
influence of the ''trade " and other prevailing winds, which give rise to slow steady 
motions which, in the absence of friction, would remain permanent even were the 
originating cause entirely to cease. The difficulty in accepting this view arises from 
the assumption that such currents would succumb to the influence of frictional 



* ' Acta Mathematical vol. 8. 

t ' Phil. Mag.,' 1875, p. 227. 

J * Hydrodynamics,' § 198. 

§ 'Physical Geography/ §§ 57^60. 



ANALYSIS TO THE DYNAMICAL THEORY OF THE TIDES. 205 

resistances in the course of a few days.* If this assumption be correct it will of 
course be necessary to invoke some more constant cause than the fickle winds in 
explanation of ocean currents, but unfortunately the causes put forward by the chief 
opponents of the wind theory, namely, the differences of density arising from 
differences of temperature, salinity, &c, though no doubt satisfying the criterion of 
being more constant in their action, seem to be equally ineffective in maintaining the 
currents against such large resistances as would be required to destroy the currents 
due to the winds in a few days. If, on the other hand, the period of subsidence of 
the free current-motions is to be reckoned rather by years, these motions could not 
fail to be excited and maintained by such causes as the winds even against the 
action of friction. 

In § 14 I have dealt with the dynamics of ocean currents on the supposition that 
they are of the nature of free steady motions (probably maintained by a variety of 
causes), and that the influences of viscosity are extremely small. A remarkable 
result is the extremely restricted character of the possible forms of steady motion as 
contrasted with the case where the ocean covers a non-rotating globe, in which latter 
case the possible forms of steady motion are to a large extent arbitrary. It is found 
that if the density of the water is uniform, the only forms of steady motion possible 
when the depth depends on the latitude alone are those in which the water always 
moves along parallels of latitude, while in general the paths of the fluid particles 
coincide with certain lines depending only on the distribution of land and water and 
on the configuration of the ocean bed. The equation by which these lines are defined 
is of an extremely simple character, and from it we could at once trace out the forms 
of the stream-lines on a chart if we had a sufficient knowledge of the configuration 
of the ocean bed. The equator will always be one of these stream-lines, and herein 
we seem to have the explanation of the fact that the ocean-currents always tend to 
set along the equator, but in other respects it is shewn that the effects of variations 
of density will seriously interfere with the simple laws which must hold so long as 
the density is uniform. 

The importance of the earth's rotation in influencing ocean-currents has long been 
recognised by physicists, but I am not aware that any previous attempt has been 
made to investigate this influence mathematically. The numerical results obtained 
in § 15 are interesting, as showing how a cause, which on a non-rotating globe could 
not give rise to any appreciable currents, may be rendered highly effective in main- 
taining currents as a consequence of the rotation of the earth. 

In attempting to account for ocean currents, the real question at issue is : How 
far are the suggested causes capable of maintaining currents against the action of 
friction ? To answer this question an investigation, either mathematical or experi- 
mental, as to the effects of friction is essential. Such an investigation I have endeavoured 

* Maury, c Physical Geography of the Sea/ § 93. 
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to supply in another paper, * but mathematical difficulties have compelled me, in 
treating of friction, to omit from consideration the important influences due to the 
rotation. We have already called attention to the fact that the free steady motions 
on a non-rotating globe are far less restricted in character than those on a rotating 
globe, while, in that the latter essentially violate what appears to be a necessary 
condition when the water is viscous, namely, that there can be no slipping at the 
bottom, it seems to me to be probable that even the limited forms of steady motion 
here dealt with would be no longer possible if the water were viscous, but that, if 
they were started by any means, they would at once give place to periodic motions 
of comparatively short period.t This conclusion has been forced on me by the 
apparent impossibility of satisfying the equations of motion of a viscous ocean on a 
rotating globe by means of slowly declining current-motions. If such should be the 
case, it follows that no stable currents can exist without variations in the density of 
the water. As however I have not as yet been able to support this view by 
anything approaching a rigorous mathematical treatment, the question must for the 
present remain open. 



§ 1. Differential Equations for the Vibration of a Rotating Mass of Liquid. 

Suppose we are dealing with the small oscillations of a mass of liquid about a state 
of steady motion consisting of a rotation as a rigid body with angular velocity co about 
a certain axis. 

Take this axis as axis of .z, and refer to a set of rectangular axes rotating about it 
with uniform angular velocity a)< Then, in the steady motion supposed, the fluid will 
have no motion relatively to these axes. 

Let u v v, iv denote the relative velocity-components at the point x, y, z due to the 
small oscillations. The actual velocity-components parallel to the instantaneous 
positions of the moving axes will then be 



u — (oy, v + <^ 3 w 



? 



and, therefore, if we suppose the amplitude of the vibrations sufficiently small to allow 
of our neglecting squares and products of the small quantities u, v, iv, the differential 
equations of motion of the liquid may be written in the form J 



* Loc, Git., ante. 

t The condition that there can be no slipping at the bottom will reduce the number of degrees of 
freedom of the system, and hence we may anticipate that certain types of motion which were possible 
before this condition was imposed will no longer exist afterwards, 

f Basset, < Hydrodynamics,' vol. I, p. 22, 



ANALYSIS TO THE DYNAMICAL THEORY OF THE TIDES. 207 



dt 



ft) (v + WX) - 0)V =r g- (Y'-p/p), 



^ + ft) (M - ft)?/) + ft)M = g~ (V - P//)), 

— - -- ( V - p/p) 

dt dz [ V VIPh 

where V denotes the potential of the bodily forces acting on the liquid, p the fluid 
pressure,, and p the density. 
If now we put 

if, = V - p/p + |ft) a (x 2 + ?/ 2 ) + const (1), 

the above equations reduce to 



ou 

-w: 2(x)V = 




& + 2ftm = 


3l/r 

"3.v 


3w 


3i/r 


a* 


9^ , 






while the incompressibility of the liquid is expressed by the additional equation 

----- -|- ~~- _j- — = , . ' . (3). 

o# dy d# 

These equations, originally given by Poincahe, # suffice, in conjunction with 
certain conditions which must hold at the boundary, for the determination of the 
four functions u, v, w 9 \jj. They are perhaps the simplest equations of which to 
make use when dealing with the oscillations of a mass of liquid of finite extent in 
three dimensions, and, for this purpose, they were first solved by Poincahe in a 
form adapted for satisfying boundary-conditions at an ellipsoidal surface, while 
additional applications have been considered by Bn\'AN,t Love^ and myself. § The 
possibility of solution in each of these cases however turns on the fact that it only 
required to satisfy boundary-conditions at a single ellipsoidal or spheroidal surface, 
whereas, in the problem presented by the tides, it is necessary to satisfy conditions 
at two surfaces, namely, the ocean bed and the free surface of the ocean. 

There is, however, a feature attached to this problem which enables us to surmount 

* ' Acta Mathematical vol. 7, p. 356. 

f ' Phil. Trans.,' 1889. 

X ' Proc. Lon. Math. Soc.,' vol. 19. 

§ 'Phil. Trans.,' 1895. 
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the difficulties arising from this cause ; the fluid which constitutes the ocean may be 
regarded as a thin layer distributed over an approximately spherical surface, a 
circumstance which enables us to reduce the number of our independent variables 
and to treat the problem as a two-dimensional one. 

Before proceeding to the transformation of our equations, let us examine them in 
the form in which they are given above. If we suppose the system is executing 
a simple harmonic vibration in period 27r/X, we may put u, v, iv, \jj each proportional 
to e ik \ and therefore replace du/dt, &c., by iku 7 &c. Thus the equations (2) give 

•\ o ^ 

ox 



iXv + 2(ou 



ikw 



3if 



dy ' 

dyfr 

dz 



Now in an important class of oscillations, viz., the tides of long period, the value 
of X will be small compared with that of o) ; while for another class of motions, viz., 
the steady ocean-currents, we must suppose X absolutely zero. In these cases, if we 
retain only the most important terms, the equations of motion take the approximate 
form 

— - 2<Wi; = ;;---, 2(OU = ;> - , = ~ -- • 

ox oy vz 

Hence, applying the operators d/dz to the first two, and making use of the third, 

we find 

du/dz = 0, dv/dz = 0, 

Likewise also from the equation of continuity, 

dio du dv 1 / cpyfr 3 3 ^ 



yjr 



dz dx cy 2(o \dx dy dy dx, 



0. 



From this we see that in the case of tides of very long period the velocity of the 
fluid paxticles is approximately the same at all points in the same line parallel to the 
polar axis, while in the case of the ocean-currents this is rigorously the case. 

Now in order to effect the transformation of the equations of motion, it has been 
assumed by Laplace and his followers that, on the analogy of " long waves" when 
there is no rotation, all fluid particles which are at one instant in a vertical line 
will remain in such a line. This assumption appears to require some modification in 
the case of our rotating system. We shall see hereafter however that the assump- 
tion in question will not lead to appreciable error, provided that the depth of the 
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water is small in comparison with the radius of the solid globe on which it resides, a 
hypothesis which will certainly be applicable in the case of the earth. 

§ 2. The Boundary -conditions. 

Before proceeding with the approximations which we propose to employ hereafter, 
let us examine the boundary-conditions to which the functions u, v, w 9 xj/ are subject 
in the general case. 

Suppose the fluid resides on the surface of a solid nucleus which is constrained to 
rotate with uniform angular velocity co about the axis of z. We introduce this con- 
straint so as to avoid the complications resulting from the reactions of the fluid 
motion on that of the nucleus. Since in the case of the earth the mass of the ocean 
is exceedingly small compared with that of the solid parts, such reactions would be 
very minute, while for most of the more important types of oscillation they would not 
exist at all In such types the problem is not affected by the introduction of the 
constraints. The boundaries of the ocean where it is in contact with the solid 
nucleus may then be regarded as fixed relatively to the moving axes Ox, Qy, Qz, and 
the condition to be satisfied at these boundaries is that there is no flow of fluid across 
them. Denoting by l 3 m, n the direction-cosines of the normal to the surface, this 
condition is expressed analytically by the equation 

\lu + 'wiv + nw ~\ = • "• • • • • • • • (4)- 

.Next, let us examine the boundary-conditions at the free surface. Let ?, m, n 
denote the direction-cosines of the normal to this surface in its undisturbed position^ 
and let £ be the distance between the displaced surface and the mean surface 
measured along the normal to the latter. Then we may equate the velocity at the 
mean surface in the direction of this normal to the rate at which £ increases ; thus at 
the undisturbed surface we have 

"hi + mv + nw ~] — 3£/d£ (5). 

Lastly, we must express the condition that the pressure at the actual free surface 
is zero (or constant). Now if dn! denote an element of the normal to the undisturbed 

surface, and p, ~-, denote the values at this surface of the pressure and its rate of 
increase along the normal, the pressure at the actual surface is 



P T d f\ f •••••••••»•• \P)i 



and this we have seen is to be equated to a constant. 

MDCCCXCVII. — A. 2 E 
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But by definition of xft we have 

■p/p = const + V + J- ct> a (x~ + 2/ a ) ra ^ 

= const + V + v' + v + i ^ (%? + ?/ 2 ) — ^ . . . . . (7), 

where V denotes the potential at x, y, z in the steady motion, v f the potential due 
to the attraction of the layer of fluid contained between the actual free surface and 
its mean position, and v the disturbing potential, which may be regarded as due to 
some external attracting system. 

Since in (6) 'dp fin is already associated with the small factor £, in calculating dp fin 
we may omit all small quantities of the order £, and thus replace this expression by 
its value in the steady motion. But from (7) we have in this case 

p/p — const + V + |o) 2 (xfi + y% 
whence 



Now. since the free surface of the ocean must be an equipotential surface, the 
resultant of gravitation, including centrifugal force, must be perpendicular to this 
surface. Denoting by g its value, we have 



<j=-*-,{r +lafi (** + !?)) 



aid therefore 






Introducing the values of p, k~, from (7), (8) into the expression (6), and equating 
the latter to a constant, we find 

[V' + v + v + |or (x 3 + 3/ 2 ) — i//] — gl — const, 

or, on equating periodic parts to zero, 

xjj = 7/ — g^ + v a . . . . * . . . * (9), 

where the bars are used to denote surf ace- values at the undisturbed free surface. 
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§ 3. Transformation of the Equation of Continuity. 

We proceed now to the transformation of our equations into a form analogous to 
that used by Laplace, dealing first with the equation of continuity. 

Let us refer to a system of orthogonal curvilinear coordinates a, /3, y, and suppose 
that the undisturbed surface of the ocean coincides with one of the surfaces y= const, 
say y = y . 

On the surface y = y take a small parallelogram PQRS, bounded by curves of the 
systems a = const, /3 = const. 



,5 



p«r>-~- 






m£ \ 




Through the sides of this parallelogram draw the surfaces a = const, /3 = const, to 
meet the inner surface of the ocean in the quadrilateral P'Q'R'S' and the distorted 
free surface in pqrs. 

The surfaces a = const, fi = const, are, of course, supposed to be in rotation with 
angular velocity co, in common with the axes Ox, Oy 9 Oz. Let IT, V, W denote the 
relative velocity-components at the point (a, /3, y) parallel to the normals to the 
surfaces of reference and in the directions in which a, /3, y respectively increase, and 
let y = y v where y x may be regarded as a function of a, /J, be the equation to the 
surface of the solid earth. 

Let h v /i 2 , ^ 3 be parameters associated with our orthogonal system of coordinates, 
such that the line-element ds is given by 

cM = dafi/hf + dp/hf + dy*[h£. 

Then the volume of liquid which flows in a unit of time across the face a = const 
of an elementary parallelopiped whose adjacent edges are Ba/h ly §/3//i 2 , $y/h s is 



U 



S/3 By 



But if in the above figure we suppose that a, /3, y are the coordinates of P, and 
that PQ = 8/J/A 3 , PS = 8a/h l9 the total flow of liquid across the face PP'Q'Q will be 



&j lii £i 
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found by integrating the above expression with respect to y between the limits y=yn 

Y = 7o- 

Thus, the rate at which liquid is entering the elementary volume PB/ across the 

face PQ' is expressed by 

fVo 
(U/A 3 A 3 ) dy. 



If in this expression we change a into a + Sa ? we shall, obtain the rate at which 
fluid is flowing in the positive direction across the face SR/ ; therefore, the rate at 
which fluid is leaving the element across the face SB/ is expressed by 



8$ \ y \V/hA) dy + 8a 80 | a { P(U/W dy 



In like manner, the rate at which fluid enters across the face PS' is 



8a (Yjh^) dy, 



Yl 



and that at which it escapes across the face QR/ is 



Sa P(V/W dy + S« SiS ^ { r°(V/A s / h ) dy 

Lastly, in virtue of the boundary-equation (4), which holds at the surface 
P'Q'JJ'S', the rate at which fluid enters over this surface is zero, while in virtue 
of (5), which holds at the surface PQRS/ the rate at which fluid escapes over the 

latter surface is expressed by 

Sec 8/3 9£ 

\ ?h dt ' 

Now the total amount of liquid contained within the elementary volume under 
consideration is constant, and therefore, if we equate to zero the sum of the rates at 
which fluid is entering over all the six faces ? we obtain the equation of continuity in 
the form 



Set S/3 



la 



Yo 



(WW d 7 



Yl 



Sa S/3 



1 

d/3 



Yo 



Yl 



V/A 3 /h) dy 



b* sp ar 

h } h z dt 



= 0. 



or 



'dt 



hK 



"3 r r ?<> 



da. 



Yl 



(U/^ 3 )^y} + ^{f"(^A)^}_ 



(10). 



So far no approximation has been made other than that involved in supposing the 
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vibrations small. If we now suppose that the depth of the water is small, so that 
To — 7i 3S a sma 'l quantity, the above equation admits of considerable simplification. 
Using square brackets to denote values at the mean surface y = y , we have, by 
Taylor's Theorem, 

u//^ 3 = [UAA] + (r-ro) ra 

whence 

(U/hA) d r = [U/^s] (yo - yi) - i (yo - jif ^(U/hh) + 



& 



(U/\h s ) 



1 • • * J 



C7o 



) • • . 



we 



Now, if y — y :i be small as supposed above, # even though ~— (U/A 2 /i 3 ) is finite 

may omit all the terms on the right except the first. This amounts to supposing 
that the horizontal velocity is sensibly uniform throughout the depth, not on account 
of the small value of its rate of variation, but on account of the small distance 
through which this variation can take effect, a supposition which is not inconsistent 
with the results of § 1. Hence, on neglecting small terms of the order (y — y^ 2 , 
we have 



and, in like manner, 



> 



Vi 



(U/h 2 h s ) dy = [U/AjjAs] (y - yi ), 



Yo 



*yi 



(Y/hJu) dy = [Y/h s h l ] (y - 7l ). 



Let h denote the depth of the ocean at the point (a, /3). Then, provided h be 
small in comparison with the radii of curvature of normal sections of the surfaces 

a = const, /3 = const, y = const, we may put 



Yo - 7i 



= ll 9 



with errors of the order of the square of the ratio of h to these radii of curvature ; 
and therefore 

f(U/M 3 )^y = A[U/A 3 ], 



yi 

fyo 

yi 



(V/Mi) dy = h [V/AJ. 



Substituting these values in (10), we find 



dt 



= - h x h 



3 






(11), 



where we have now used bars to denote surface values. 



* The standard of comparison is considered in the next section. 
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If we suppose that the free surface of the ocean is a spheroid of revolution about 
the axis Oz, it will be convenient to refer to a system of spheroidal coordinates /x, </>, v 
i elated to x, y, z by the equations 

x — c ^/(l + v 2 ) y/(1 — jl 2 ) cos <£, 
y = c ->/(! + *> 2 ) .v/(l — /x, a ) sin <£, 



The line-element ds for this system of coordinates is given by 

i. ~T" V" 



1 — /^' 



whence, if we identify fi } cj>, v with a, /3, y respectively ? we have 



1 _ C^(lfi + fi*) 



1 



= c v/(l + **) v/(l - //) ; 






and, supposing that v = p is the equation to the free surface, the equation (11) 
becomes 



3/5' 



1 



a 



c vv 3 + v) d ^ 



. ^ (1 IJ?) h U J « x/(v + 1) 3<M 7(1 - /**) 



We have already neglected on the right small terms of the order h compared with 
those retained ; we now propose to make the further hypothesis that the spheroidal 
surface of the ocean is of small ellipticity e. In this case c will be small and v 
large, in such a manner however that cv Q is finite and equal to the polar radius a ; 
further l/z/ 3 will be approximately equal to 2e. Hence we find 



dt 



; [I k<> -<•■>*■*}+£{•<£ 



H?) 



'.« • » 



(12), 



where the terms omitted on the right are of order h and of order e compared with 
those retained, 

S 4. Transformation of the Dynamical Equations, 

Let # denote the inclination to the axis of z of the normal to the surface v ™ const, 
through any point ; then the direction-cosines of this normal will be 



sin 6 cos fa sin sin ^, cos $ 
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and the direction-cosines of the normals to the surfaces /x = const, <f> = const, will be 

— cos 6 cos <)>, — cos sin <f), sin 6 9 
— sin </) 9 cos <f>, 0, 

respectively. Hence we have 

U = — (u cos $ ~\- v sin (/>) cos -{- w sin 0, 

V = ?; cos (f) — u sin <£, 

W = (^ cos </> + ^ sin <i) sin # + w cos #, 

from which we obtain 

— - U cos + W sin = z^ cos <£ + i? sin <£. 
Again 

/^ — = — =-*- cos 9 cos t/ — -— sin 9 cos # + ^7 sm #> 
and therefore from (2), we find 

h y ~- = — / --x- — 2covj cos <£ cos # — f ^7- + 2am ) sin <j> cos # + ^7 s * n $ 

3U . rt TT » 
= -57 + 2cov cos 6>. 



Similarly 



ai: 



nd 



/ley zr~ = '- COS — -^~- Sill 

2 30 3y ^ 3# r 

= ( -^7 + 2gw J cos — ( jr 2<wy } sin 

= ^ + 2co ( W sin - U cos 6), 



/i-o -~- = ~- sin # cos + -7- sin o sin + -~- cos u 
dv dx T dy r dz 



2o)v\ sin # cos <£ + ( ^7 + 2am] sin sin -f ^7 cos 



/ die 
37 " 

77 — 2cov sm 0. 

ot 



Allowing for the differences in the notation, the three equations just obtained agree 
with those given by Professor Lamb/''' If we suppose U, V 5 W, \jj, each proportional 
to e iK \ they may be written 



# i 



Hydrodynamics,' p. 344. 



210 



MR. S. S. HOUGH ON THE APPLICATION OF HARMONIC 



i\JJ +.2e>V'cos0 



h — -*■ ' 



iXV - 2o) (U cos 0-W sin 0) = h 



2 



i\W-2o>Vsin0 






7, <tk 

I I/O r\ « 

"dp 



> 



. (13). 



The equations in the form we have just written them will hold good whatever be 
the depth of the ocean or the ellipticity of its' surface. We now proceed to introduce 
approximations similar to those of the last section. 

In the first place we have, as in § 2, W == p- , and this by (12) we see is of the order 

h/a compared with U or V. Hence ? omitting terms of the order h/a, and of order e 3 
compared with those retained, the equations (13) take at the surface the approximate 
forms 



iXU + 2o)V cos 6 



iXV — 2ft>U cos 6 



2o>V sin 6 



\/ (1 - H?) If 

O/JL 






a 



1 



df \. 



a sj (1 — /i 3 ) .90 

9-^ 
9^?/ ' 



r 



• ( 14 ) 



where, in conformity with the notation of \ 2, we have denoted by dxjj/dn the rate of 
variation of iff in the direction of the normal to the surface of the ocean. 

From the equations (14) it appears that d\Jj/dn' is a quantity of the same order of 

. , ■ - : 9 

magnitude as \jj/a; also if we apply the operator "k^r- to each of the equations (13), 

we shall obtain equations which enable us to express dXJ/dn\ 9V/3ft', dW/dn in terms 
of U, V, W, and the surface-values of xj; and its differential coefficients. A little 
consideration will show that in general 9U/3??', 9V/9?i' 5 dW/dn must be of the 
order U/a, V/a. 

Now the approximations introduced in the last section will hold good provided 
that we may neglect 

r 9 i r 9 

(yo - yi) I g^ ( U /V 3 )J and (y - Yl ) g- (V//iA) 



in comparison with [Q//i 2 ^ 3 ], [V/Agfrj 
neglect 

8 



; or, in our present notation, that we may 



a 



h ^ {cV V + ^) V (1 - /* 3 ) U) and A ^ 



C 3 <> 3 + ^8) V 



dn' [ v /(l_^) v /(l + ^ ) 
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o 



17 



in comparison with 



TJeV ("o 8 + i*')v'(l -/*'). V 



c 3 (d 2 + yU, 3 ) 



^(l-^V(l + "o ! ) 



But we have seen that dU/dn, dV/dn' are of the order "U/a, while 

3»' i c v V" ^ ** '' 9 a ' 



v/ (1 + ^ 



are of the order tfvja. 
Hence 



a 



a»' 



c 3 (y 3 + fj?) V 



lV /(l-^ )%/(1+ ^ 



are both of the order 



U 



r * o 



a 



O , 



The approximations will therefore be admissible, provided h/a is a small quantity, 
that is, provided that the depth of the ocean is small in comparison with the radius 
of the solid earth, a hypothesis as to the validity of which there can be no doubt. 

Returning now to equations (14), and solving for II, V we find 



U(X» - 4«* cos*0) - - iX ^^) ^ + _*2™1_ d ± 9 

a dfi a yj (1 — fj?) dcj) 

V (X 8 ~ 4c/ cos* 0) = - 2« cos ^ (1 " ^ |* - - -7#~ 9 - 3 # 



But we have rigorously 



cos * 



/* 



a/ (^o 3 + 1) 

VW+1*) 



and therefore, with errors of the order of the ellipticity, we may replace cos 6 by /*. 
Hence, finally, we obtain the values of U, V with errors of the order h/a, c compared 
with their true values in the form 



U = 



i\ *J (1 — fj?) d\]r 



m m mum 



9, 



(Vfjb 



d-ifr 



-\ 



a (\a - 4© V) 3/* a \/ (1 - /^ 2 ) O* - 4ftJ V 3</> 



V 



MDCCCXCVIL— A, 



2&)/^ y/ (1 — ^3) 3-^r 



^A, 



9-v/r 



• • 



iv/d-^)^- 4»V) 9^> ! 



a 



2 * 



(15). 
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Substituting these values in the right-hand member of (12), we obtain 



ika 2 £ 



3 f i\h (1 - ij?) 3^ 



2a)/jih 3t|t 



dfM 



9.. 9, 



X 2 — ia^fM 



d/J, X 2 — 4&) 3 /^ 3 dcj) 



+ 



2 60 {lib d^fr 



+ 



iX/fc 



3\|f 



3</> [X 2 - 4®V 2 3/a ' (1 - /^ 2 ) (x 2 - 4ft>V) 8</> 



e a 



(16). 



Hence, provided that X be not equal to zero, we have 



a^ 



3 f h(l — fJ?) 3^ 



2eo/jbh 






d/i\\* - 4o>> 3 3/4 iX (X 3 - 4«V) 3<£ 



+ 



0(f) 



2cojjl7i 



v^jr 



+ 



h 



difr 



iX (X 2 — iw^/A?) dfX, (1 — yt6 2 ) (X 2 — 4:6) 2 yU, 2 ) ?<^> 



• • (in 



This equation, in conjunction with the pressure equation (9) of § 2, serves to determine 
*/>, £ in terms of /x, <£. It is equivalent to the well-known equation used by Laplace 
in the ' Mecanique Celeste.'* Omitting from consideration for the present the types 
of motion defined by X = 0, we propose in the present paper to discuss only those 
solutions which are symmetrical with respect to the axis of rotation. In order that 
such solutions may exist, we must suppose that h, xfi are independent of <f>; the 
equation (17) will then reduce to 



3 \ h(l-~ ^ ) df] 
dfi 1 p - fi 2 3/i 



4a 3 tt> 3 £ 



• 



(18), 



w 7 here for brevity we have put \/2<o = f. 



§ 5. Integration by Means of Zonal Harmonics. 
Suppose that £ is expressible as a series of zonal harmonics of the form 



n =oo 



£ — s c«p w (/x)« 



n — \ 



Neglecting the ellipticity of the surface, we may at once write down the value at 
the surface of the potential due to this distribution ; we have, namely, 



n=\2n -f- 1 Vl ' 



* Part I., Book IV., § 3, 
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where the density p is expressed in gravitational units, 
density of the earth as a whole, including the ocean, 
degree of approximation, 



But if & denote the mean 
we have, with the same 



and therefore 



W=QO 



v' = X 



g = |- 7rcr«, 



S P 



n—l 



(2n+l)<r gGnPn ^ 



As this only involves the ratio p/cr, it is independent of the unit of mass employed. 

Next suppose that the surface-value of the disturbing potential can be expressed 
by means of the series 

2y*P« (/*)• 

Then equation (9) of § 2 gives 



^ — v' — gi^v 



•6p 



g ( l _ ::r ) Q 

J l (2n + 1) a- 



n 



7 



n 



or, if we write for brevity, 



9* = 9(1 — 



3 P 



r = 



n 



we have 



Now, from the equation 

d 

cI/ju 



dP t 



V n (fi) 9 



(2n+l)aJ 9 ^ 9 . 



i/r = 2l\,P« (/x) . % 



• 9 O 



» • * • • 



(l-^)^j +»(n+l)P. = 



(19), 



(20). 



which defines the zonal harmonics, we find 



7? (n + 1 ) P» dft. 



Hence, if we integrate (18) with respect to ft, we obtain 






AaWSC n ^F n dfi 



a 



rfp« 



= A - 4a V (1 - ^) t -r^-r. - r n , 

\ r 7 7t (ft + 1) dyU, 



where A is an arbitrary constant, which may be seen to be zero by putting //, = ± 1. 
Therefore 

2 F 2 



09f) 

i~i £j \J 
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h (i - ^) y= ~ AcM i 1 - 1^) S 



C 



» 



rfP 



«■ (« + l) 



</ 2 -" ! »t ■ ■ < 31 >. 



But. by well known properties of the zonal harmonics, we have 



a -/••>£ = 



— n(n + 1)1 P«cfyx, 

_ ft (ft + 1) p /dPn+i _ ^~~1 

2% + 1 J \ rf/A rf/A / 

^ (w + 1) , p p V 



) d/x 



271 + 



no arbitrary constant being necessary since both sides vanish when jx '== 1 ; and 
therefore 



(/ 2 - /* 2 ) 






(/* - 1) 



(7 3 -i) 



dP„ %(% + 1) 

rf/A 2% + 1 

dP n n(n -\- 1) 

cZ/z 2^ + 1 



(l»+l *— P»-]) 



riP 



^ 



dP« rfP 



w 



whence 



«-2 



2w 4- 3 \ ^ clfjb J 2% —l\dfi dfi j ^ 



(P - ^) 



d? 



n. 



n (n 4- 1) dl 



»+2 



tf^t 



(271 + 1) (2n + 3) djLL 



+ /» - 1 4 



2ri(n 4- 1) 



dP 



n (n 4-1) 



tZP 



w— 2 



(2n - 1) (2w + 3)/ ^ (2w - 1) (2n +1) rf^ 



This relation will hold good when n = 1, provided we replace dP^/cfyx by zero. 
Thus the right-hand member of (21) is equal to 



r»=co 



r 



4a 3 <» s (l— jt 2 ) X 

= 1 



C«_ 



«-2 



_(2rc-3)(2?i-l) 



^« 



n (n + 1) "*" (27i- 1) (2w 4- 3)/ "^ (2% + 3) (2ra 4- o) J 7/7 



The left-hand member may in virtue of (20) be written in the form 



1 






Equating the two members and dividing by 4co%^(l — /x 2 ) ? we obtain 



4coW* * dfi 



d? 



c tt _ 



dp l(2w-r-3)(2w-l) 



C, ( -•?— i + 



2 



T" 



^M 



+2 



n(n 4- 1) ! (2w - 1) (2n + 3)/ ' (2n 4- 3) (2ro 4- 5) J 9 
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Hence, provided h be constant, the two members will be identical, if for all positive 



integral values of n we have 



C 



n—9. 



(2ft - 3) (2ft - 1) 



CA^'A + T, 








^n + % 



n(n + l) * (2rc-l)(2rc-3)/ ' (2rc + 3)(2w + 5) 4a> 2 a 



= aZ^s ( 22 )> 



it being understood that O = 0, C_, = 0. 

This is on the hypothesis that the depth is constant ; a more general hypothesis 
would be to suppose that h is of the form 

k +• 1(1 — jLt 2 ) 

where k and I are constants 

Assuming this form for h, the left-hand member of (21) becomes 

rlV rlV 

h(\- it) tr n ~ + 1 (l - ft) sr„ (i - ^) " 



d/A 



which, by the properties proved above, is equal to 



d 



^ 



7, 

It 



(1 -ja 2 )sr w 



dp 



w 



dfju 



— 1(1 — /x 3 ) 



^ s 



_0L-2)_(rc-l)_ 
_(2rc- 3)0^-^1) 



■p ^^ 



_Jn («+ 1) _ (ft + 2) (ft + 3) 

(2ft - 1) (2ft"+~3) " + (2ft + 3) (2m + 5) * 






dP. 



Identifying this with the right-hand member, we obtain 



C 



M-2 



(2ft - 3) (2ra - 1) 



o 



4» 2 a 2 4w 2 a 3 



/ 3 -l 

71 (71 + 1) 

i\(2%-3)"(2%-i) 



"T" /k 



(2w - 1) ^ + 3), 
2^0+1) 










n+2 



(2a, + 3) (2n + 5) 



v J- 11 O 



n-% 



X * "+" (2^ + 3)(2^H-5) lf ^ 2 i \^ A /' 



On introducing the values of r„ from (19), equations (22), (22a) may be written 



W 



(2?i - 3) (27?, 



1) "" n \n (n + 1) + (2w ~T)(2^ +1) 









C 



rc+2 



(2w + 3) (2w + 5) ~" 4<»%< 



• a o t. e » 



• • (23), 



C 






&Y» 



"(271 -3) (27?, -1)~"~~ 






"T" O. 



(/« - 1) . 2 {1 - ft (ft + 1) lg»l±afia*} kg„ 



n (ft +1) (2ft - 1) (2ft + 3) 

1 - (ft + 2) (ft + 3) Igml^aW 



I 



n+i [ (2ft + 3) (2ft + 5) 

(ft-2)(ft-l) 



4 W % 3 



4o) 2 a 3 4o)% a 1. (2ft- 3) (2ft- 1) 



; y»-» 



2ft (ft +1) __ , (ft + 2)(ft+3) 

on y» "i" /o„ i m\ /•>„ ,k\ y»+2 



(2ft-l)(2ft+3) / " ' (2ft+3)(2ft+5) 



• ( 2 3a) 
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The law of variable depth which we have assumed appears to be the most general 
law which will lead to a difference-relation connecting the successive (7s of order not 
higher than the second. We shall confine ourselves chiefly to the case where the 
depth is uniform, but the following remark with respect to the more general case 
seems worthy of attention. 

If we put for brevity 

c. __ 1 — n(n + 1)^/4g)% 3 _ 1 ~" n ( n + 1) Igjil &aPa? 






n 



and suppose all the ys zero, so that there is no disturbing force, equations (23a) 
mav be written 

- LA + 17 A — 0, - LA + %0 4 = 0. 1 

f A — L3C3 + ??A = 0, f A — L 4 C 4 + %C 6 = 0, 1 ^ ^ / 24 y 

f 3 C 3 — L 5 C 5 + t? 5 C 7 = 0, f A — LA + %A = 0, 



^ 



Now, suppose Z in the expression & + ?(1 ■— f^ 2 ) f° r the depth is of the form 
where r is an integer, and for greater definiteness let us suppose that r 



v (r + 1) g r 
is even. 



Then ir = and rj r „ 2 = 0, and therefore the equations (24) will all be satisfied if 

— L 2 C 2 + 17 A = 0, 

?A - L A + 1? A = 0, 



* *. s # a b 



f1._2Cv.-2 — LA + ^rC r+2 == 

— L r + 2 C r+ 2 + ^M-2U-f4 ^ 
§V+2^r + 2 "™" -* J r + 4^V+4 J ^r+4^ + 6 ^ " 



3 



I 



? 



and all the C's with odd suffixes vanish. 

Further, these will be satisfied if C v+2j O r+4 , . . . are all gero, provided X is a root ol 
the equation 



ANALYSIS TO THE DYNAMICAL THEOEY OP THE TIDES. 



223 



L 2 , 7) 2 , 0, 0, ........ . 



o, 



£ 4 , L*6> ^o 



L r _ 6 , 


Tjr-Qi 





br-65 


— L. r __ 4 , 


^r-4 


o, 


€V-4? 


-L r _2 



0. 



It follows that there exist certain types of free oscillation for which all the values 
of C with suffixes greater than r are zero. For these types the height of the surface- 
waves will be expressible by a finite series of terms terminating with a term 
involving P r . 

In like manner, if the disturbing force be derivable from a potential function of the 
form of a second order harmonic, the equations which determine the forced 
oscillations are 



L ?A + %c 4< = Tztk + 



4 ly 2 



4oA^ 



7 4a)%3 ' 



£A — L 4 C 4 + tjJJq 



2 . 3 lyo 
577 4«% 3 ' 



£A — k 6 C 6 + >?6^8 = °> 






+ 2 



= 0, 



-L | r4-2^'/'+2 'F 7 ?r+2^r + 4 — 0. 



# » 



• -•* 



If we suppose C r+2 , C r+4 , . . . all zero, the first r/2 of these equations will serve to 
determine C 2 , C 4 , . . . C> in terms of y 2 , while the remaining equations will be 
satisfied identically. Thus the forced tides for the law of depth in question will be 
expressible by a finite series of terms terminating with a term involving P r . 

This general law does not hold when r = 2, owing to the presence of a term in the 
right-hand member of the second equation. 

The fact that for these laws of depth the tide-heights could be expressed by finite 
series, instead of by the infinite series usually required, was originally proved by 
Laplace in the ' Meeanique Celeste/''' 



* Part I., Book IV., § 5, 
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§ 6. The Period-Equation f 07* the Free Oscillations. 

Returning to the case of uniform depth, we see that the equations (23) divide 
themselves into two groups, in one of which only even suffixes and in the other of 
which only odd suffixes are involved* We therefore conclude that the types of 
oscillation divide themselves into two classes, in the former of which the height of the 
surface- waves will be expressible entirely by harmonics of even order, and in the 
latter by harmonics of odd order alone. An exactly similar treatment is applicable 
to each of these classes ; we shall therefore select for discussion the former set 5 
ontenting ourselves as regards the latter with merely stating results. 

Denote by L* the expression 



c 



/»-l 






9, 



n (n + 1) (2n — 1) (2n + 3) 






(25). 



Then, putting all the y's equal to zero, the types of free oscillation will be deter- 
mined by the equations 

a 



-. 



^2^2 T" hr n — 



5.7 



7.9 



Oa-Lj, -j- 




11.13 



= 






(26). 



C 



«-2 



(2% - 3) (2n - 1) 



G fl L H + 






(2n + 3) {2n + 5) 



— 



_y 



At first sight it might appear that whatever be the value of X these equations will 
serve to determine C 4 , C 6 , . . . in succession in terms of C 33 whereas we know that this 
should only be possible for certain determinate values of A corresponding to the 
different periods of free oscillation. The manner in which these values of X are to be 
determined involves arguments similar to those used by Kelvin* in justification of 
the procedure of Laplace with reference to the forced oscillations after it had been 
attacked by Amyt and FerrelJ. 

From equations (26), we obtain by actual solution 

a 



7.9 



VVoJ-Jc 



C, 



l . J. llilO 



Z & 



Go 



iJ.T, 



1 



. / 



1 

7.~9 

Li.4. 



* ' Phil. Mag.', 1875. Gf. also an analogous problem treated by Xiven, 6 Phil. Trans.,' 1880, Part I, f 
p. 133 et seq. 

f " Tides and Waves," § 3. 

% " Tidal Researches " (U. S. Coast Survey 1873), 
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■Of 



and in general 



(_l)«/2 



^»f 2 



7.9.11 .. . (2» + 5) 



c 



3 



L 



2' 



7.9' 



0, 



1 



5.7' 



0, 



La, 



1 
11.1 o 



4,) 11 1 «>3 



d 



o 



i 

9.11' 



t ^ 

1J6 > 15. 



17 



0, 0, 



o, o, 



-L 



0, 



n - 2 ' (2»-l)(2» + l) 



(2» - 3)(2» - 1)' 



or, denoting the determinant which multiplies C 2 in the last equation written down 

by A„, 

C s + 2 =(-l) M/2 7.9.11 . . . (2»+5)A H C 2 , (> = 2,4, 6, . . . ) 

Now the determinant A„ is an algebraic polynomial of degree n/2 in f 2 . - If therefore 
we equate it to zero, we should obtain an algebraic equation of degree n/2. If f z 
has as its value any of the n/2 roots of this equation the first n/2 of equations (26) 
will be consistent, while C H+2 will vanish. By increasing the value of n we shall 
approximate more and more closely to the case where an infinite number of such 
equations are satisfied, while we shall impose an additional condition on the C's, 
viz. :— that at some stage one of them must vanish. Though in general in the actual 
motion none of the quantities C 2 , C 4 , . . . are zero, they are however subject 
to an important restriction, namely, that the series C 2 , C 4 . . . must form a 
converging series, and therefore we must satisfy the equation 



V, ~ CO 



The latter equation may be regarded as the period -equation for the free oscillations. 
It follows that as n is increased the roots of the equation A„ = 0, which make 
C„ + 3 == 0, must approach closer and closer to certain definite limiting values, which 
correspond to the different periods of free oscillation, and that the series C 2 , C 4 . . . 
calculated in succession from equations (26) can only form a convergent series when 

/has one of these values. 

Now we see that A n = is the equation obtained by eliminating the C's from the 

pet of equations 
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C 



~~ C 5 L 3 + y~q ~ °' 






[2r-3)(2r- 1) 



C,L, + 



C 



r+2 



(2r + 3) (2r + 5) 



• • a • » • e # 



0«_ 



»-2 



(2ra - 3) (2% - 1) 



C W L W 



0, y . . . . (27) 







From these we obtain 



£> - g/y r _ 



(2r - 3) (2r - 1) 



Lr 



(2r + 3) (2r + 5) 



WW+ 2 



• XJ r "— • 



and therefore by successive applications 



C r ^/O r =(2r-^3)(2r- 1) L, 



(2r+ l)(2r + 3) 3 (2r + 5) 



C/./C r + 



(2r + 1) (2r + 3) 



(2r + l)(2r + -3) 8 (2r +5) 



M- + 2 
1 



(2rc- 3) (2ft- l) 3 (2;*,-f 1) 



L M 



Thus the eliminant of equations (27) can be expressed in the form 



t _ IiL _ ^112 9.11 3 .13 

^ J 2 n in V t '""" 



e e e 



(2ft - 3) (2n - If (2ft 4- 1) 

■■■■■■>■>■■•■■■■"■■■■■■■■■■■■■■■■■■■ 



We therefore see that the roots of the equation A„ == are the roots of 



5.7 3 .9 9. II 3 . 13 



(2ft - 3)(2ft--l) 3 (2ft-f 1) 



'3 



L 



0. 



6 



hi 



This form for the equation A u = has an advantage over the determinantal form, in 

that it enables us at once to proceed to the limit when n is made infinitely great, 

and thus to express the period-equation for the free oscillations by means of the 

transcendental equation 

1 1 



L *~ L 



5.7 3 .9 9. II 3 . 13 



'4 



L 6 ~ . . t ad inf. 



~ U • » . . . {JiOj. 



ANALYSIS TO THE DYNAMICAL THEORY OF THE TIDES. 227 

We may however obtain a number of alternative forms for our period-equation. 

From (27) we find 

1 

C f ± a/fi. T _ (2r-3)(2r-l) g (2r + 1) 

(2r + 3) (2r + 5) ~~ ' q/C > ._ 2 

(2r - I) (2r + 1) 

and therefore by successive applications 



Cv + s /0, _ T (2r - 3) (2r - l) 2 (2r + 1) 5 .7^9 



(2r + 3) (2r + 5) ' L,._ 2 - ... - L 



2 



Thus, we have 

1 



(2r - 3) (2r - iy (2r + 1) (2r + 1) (2r + 3) 2 (2r 4- 5) 

q/Qf-g _ r /c r+2 

(2r - 1) (2r + 1) (2r + 1) (2r + 3) 

1 11 



(2r - 3) (2r - l) 2 (2r + 1) (2r - 7) (2r - 5) 2 (2r - 3) 5 . 7 2 ■ 9 



(2r+l)(2r+3) 2 (2r+5) (2r+5)(2r+7) a (2r + 9) (2n-3)(2n-7f(2n+l) 

1 1 t t r * 

This is an alternative form for the equation A* = ; by making n infinite, we 
obtain as an alternative form of the period-equation 



T _ (— ~~ 3 ) ( 2r Z 1 ) 2 ( 2r ± *) (2r - 7) (2r - 5) 2 (27- - 3) 5.7 3 .9 

Jm ~ J t y ~mr *^- -y 

1 1 



(2r + 1) (2r + 3) 2 (2r + 5) (2r + 5) (2r + 7) 3 (2r + 9) , . 

where r is any even integer, 

§ 7. Numerical Solution of the Period- Equation. 

The method I have used to solve the above equation will perhaps be best explained 
by giving a numerical example in detail. 

Taking for the ratio of the mean density of the earth to that of water the value 
given by Boys, # we deduce 

* * Roy. Soc. Proc,,' vol. 56, 1894, p. 132. 

2 a 2 
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p/o- — -18093, 



whence 



9m - 


•89144 


9 12/ 9 — 


•97829 


9 J 9 = 


•93969 


9n./9 = 


•98128 


9J9 = 


•95825 


9^9 = 


•98355 


9b/9 = 


•96807 


9i?l9 ~ 


•98533 


9io!9 -- 


•97415 


9J9 = 


•98676 



the values of g n jg approximating closer and closer to unity as n increases. 

Next take hg/ico^a? = -^g-, which corresponds to a depth of %-£$■ q of the earth s 
radius, or about 7260 feet. With this value of the depth we find 



L, 



I* 



L„ = 



1 



9 9 



4« 3 

/X 3 



4.5 \4ft) 3 



1 



X 3 



L 



8 



-Li in 



6.7 \4o> 2 

1 /jv 3 

879 \46) 2 

1 /X 3 



1 + '07295, 



— I + '002482, 



1 - '011835, 



1 



•017184, 



'10 



10.11 \4e» s 



1 — -019777, 



-Lif 



1 






12 



12.13 \4 



1 



•021237, 



or 



Ul/I - — 



X s 



'14 



L 



14 . 15 \4« 3 
1 /\ 3 



16 



T _ „ 



l . 17 \4 W 2 
1 /\ 3 



18 18.19\4w 3 



L 



l 



\ 3 



20 



20 . 21 \4ffl 3 



1 — -022143, 



1 - -022746, 



- 1 — -023168, 



•023476. 



Introduce for brevity the notation 



H = 



K 





1 


(2n H 


h 1) (2m + 3) 3 (2m + 5) 


1 


(2m - 


- 3) (2m, - l) 3 (2m + 1) 



1 



(2n 


— 


3) (2m - 


1)' 


(2m 


+ 


1) 






1 


2 








(2n 


+ 


1) (2m + 


3) 3 


(2m 


+ 


5) 



B « » 



n 



L 



» 



■k»+2 



1 

JL/n 



arf Wl/. 



The period-equation (29) may then be written 



JL-Jfy """" J L 



#-2 



jl\_^ 1 o * yj. .».'.* ■.«.'■•» • • • • 1 Ov /• 



Suppose now that A, 2 -/-4&>? has a value found by equating to zero one of the quantities 
L, say for example L 8 ; putting L 8 = 0, we obtain with the numerical values given, 
above ' ' - ; - : : 






4co 



,3 



= 2*23726, 



L 10 •— 


- -008529, 


L i2 = 


- -013306, 


L 14 = 


- -016251. 
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and with this value for X 3 /4o) s we find 

L = -27916, 
L 4 = -064345, 
L 6 =°'017624, 

The value of log — -7~7r7o~T ~^\* Vo — 7"T\ * s 

° (2n + 1) {2n + oy {2% + 5) 

for n= 2, 4-6566,. for n = 8,6*8898, 

„ n =.4, 5-8490, „. n = 10, 6*5564, 

„ n=6, 5-3034, „ n = 12, 6'2769. 

With these values we find for the successive convergents to the continued frac- 
tion H 6 

•001141, -001217, '001219, 

while the successive convergents to the continued fraction K 10 are 

— "•000910, - '000939, - -000940, 

It will be observed that these continued fractions converge with great rapidity ; so 
long as the depth of the ocean is not less than that we are here using, I find that 
when X 2 /^ 2 has a value in the neighbourhood of a root of the equation L„ = 0, the 
continued fractions H w _ 3> K„ +2 are represented without sensible error by their fourth 
convergents, while in many cases the second convergents will form a sufficiently 
accurate approximation to their values ; this rapid convergence of course greatly 
facilitates the numerical computation. In practice, the simplest method of evaluating 
the continued fractions is to assume that, for a sufficiently large value of ■% K„ = 0, 
and then to compute K*_ 3 , K^ i5 &c, in succession from the formula 



v j2n-l) (2n ±- 5)* (2n - 3) # 



Thus, in the present instance, we may put K 16 ■ = 0, and deduce 

I°g--Ki4 = ^4-06, log K 12 = 7*4*436, log K 10 = h4*9730, 

whence, as above, 

K ]0 = - -000940. 
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In like manner, in order to evaluate H w _ 2 , we assume that, for some sufficiently 
small value of r, H r = 0, and then compute H r+2 , H f+4j , • • • * H*_ 2 in turn by 
means of the formula 

Lf-t-2 "~" i-ij- 

In the present case we find 

log H 2 = 3*2107, log H 4 = 3-0516, log H 6 = 3*0860, 
whence 

and 

L 8 - H 6 - K 10 = - '001219 + '000940 

= — '000279. 

This being a small quantity we conclude that there is a root of the period-equation 
differing but slightly from the value assumed for X s /4co 3 , namely, 2*23726. A closer 
approximation will be found by using this value in H 6 , K 10 , and again equating 
L 8 — H 6 — K 10 to zero ; in other words, by putting 

L 8 = + '000279. 

The second approximation to the root is therefore given by 

X 2 /4g>* = 2-25735. 

Taking this value, and proceeding as before, we find 

L 8 - H 6 - K 10 = '000279 - -001183 + '000961 

= -000057. 

We have now found that, when X 2 /4o> 2 = 2*23726, 

L 8 -H 6 -K 10 = - -000279, 
and when X*/4g) 3 = 2*25735, 

L 8 - H 6 - K 10 = + -000057, 
whence, by interpolation, we conclude that 

L 8 — H 6 jv 10 ^ 0, 

when 

X s /4w* = 2'25394. 

In general we shall at this stage obtain a sufficiently close approximation to the 
root sought, as may be verified by actual substitution. Should however great 
accuracy be desired, we may re-start the computation, using the value already found 
as a first approximation, and so continue until the desired degree of accuracy is 
attained. The number of cases in which I have found a repetition of the process 
necessary is however extremely limited. 
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By the method here sketched I have calculated the first six roots of the period- 
equation for four different depths of the ocean corresponding to the values 
*¥o> lufj *To> ih respectively, for hgjiofia 1 . These depths are equivalent to about 
7260, 14520, 29040, 58080 feet respectively, and the results are embodied in the 
following; tables : — 

X iVoLE J.. 



n 

n 
n 
n 
n 
n 



2 
4 
6 
8 

10 
12 



n = 


2 


n = 


4 


n = 


6 


»i = 


8 


w = 


10 


w = 


12 



n 


— 


2 


n 


— 


4 


n 




6 


n 


— 


8 


n 




10 


n 


— 


12 



n 
n 
n 

n 
n 

n 



6 
8 

10 
12 



A pproximate value 

of X 2 /4w 9 , computed 

from equation 

L n = 0. 



Corrected value 

of\ 2 /4oA . 



Period of 

oscillation expressed 

in sidereal time. 



h. 



m. 



Corresponding 

period wben there 

is no rotation. 



Depth 7,260 feet (hg/laPa* = -£>) ; p/tr = '18093. 



•56230 
•95036 
1-4971 
22373 
3*1755 
4-3130 



♦44155 
•96357 
1-5224 
2-2539 
3-1867 
4-3209 



18 


3-5 


12 


13-5 


9 


43-5 


tj 


59-6 


6 


43-3 


5 


46-4 



Depth 14,520 feet (hg/4a>*a* = ■&) I pl<* — '18093, 



•69600 
1-4202 
2-5032 
3-9798 
5-8544 
8-1283 



•62473 
1-4368 

2-5168 
3-9882 
5-8600 
8-1322 



15 


11-0 


10 


0-7 




33-8 


6 


0-5 


4 


57-4 


4 


12-5 



Depth 29,040 feet (%/4o> 2 a 2 = -f -); p/cr = -18093, 



•96344 

2-3599 

4-5155 

7-4649 

11-2123 

15-7588 



•92506 

2-3707 

4-5224 

7-4691 

11-2150 

15-7609 



12 


28-6 


ti 


47-6 


5 


38-6 


4 


23-5 


3 


35-0 


3 


1-4 



Depth 58,080 feet (%/4<o 2 a 2 = \) ; p/cr = '18093. 



1-4983 

4-2393 

8-5402 

14-4352 

21-9275 

31-0202 



1 -4785 

4-2453 

8-5437 

14-4371 

21-9293 

31-0212 



9 


521 


5 


TiJ To 


4 


6-3 


3 


9-5 


2 


33-8 


2 


9-3 



h. 



m. 



32 


49 


17 


30 


11 


58 


9 


5 


7 


20 


6 


9 



im*J 


12 


12 


23 


8 


28 


6 


26 


5 


11 


4 


21 



16 


25 


8 


45 


5 


59 


4 


33 


3 


40 


3 


4 



11 


35 


6 


11 


4 


14 


3 


13 


2 


36 


2 


10 
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By a comparison of the 2nd and 3rd columns it will be noticed that in most cases 
the roots of the frequency-equation are given at once with a fair degree of accuracy 
by simply solving the equations 1% = 0, and that this approximation improves the 
greater n becomes. In the fifth column I have given the periods of oscillation for an 
ocean of the same depth when the rotation is annulled, calculated by means of the 
formula 

X 2 n (n + 1) hg n 

where co now denotes a constant such that rr/co = 12 hours. It will be seen that the 
approximation obtained by omitting the rotation continually improves with increasing 
values of n 9 but in no case will it lead to as accurate a result as the formula 

1 = n (n +1) 



4o> 2 v i / (_4or% 3 (2n - 1) (2n + 3) 

For instance, taking the case hgj&uPa* = -^o, w, = 8, the error introduced by using 
the first formula amounts to about 14 per cent* 5 whereas the second form gives the 
frequency with an error less than one per cent, of its true value. 



- § 8. Unsymmetrical Types. : 

An exactly similar method of treatment is applicable to the types which are 
represented by a series of harmonics of odd order ; the period-equation for these types 
is given by 

T TT — TT — 



where n now denotes an odd integer and H w , K n denote respectively the continued 
fractions 

1 11 



( 2n + 1) (2n + 3) 2 (2n + 5) (2k - 3) (2n z I) 2 (2rc + 1) 3 . 5 2 . 7 

1 1 



(2u - 3) (2% - l) 3 (2» + 1) (2« + 1) (2w + 3) g (2w + 5) 

n - ^ 1 — — T Ttrn-r— ni w wiB ■ nm nini n ■■■■ imiinimnwi w — i|wwmiwiiiiwi«wrmiiiiiiiiiiiiwiiMiniiin a mca^ammaeammm&&2mammm& mmsttammBm»am 

Treating this case in the same manner as the last, I have found the first six roots 
and the corresponding periods of oscillation for the four depths employed as follows :— 
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X A.I31jE JLl* 






! 


Approximate value 

of \ 2 l4w\ 


Corrected value 
of \ 2 /4w*. 


Period of oscillation. 


Corresponding 

period 

without rotation. 








h. m. 


h. m. 




Depth, 7,260 feet {hgjio^a? = - 4 \) ; p/cr = 18093 


• 


71 = 1 


•24095 


•15491 


30 29-3 


59 17 


n= 3 


•74340 


•70890 


14 15-2 


22 49 


n = 5 


1-2002 


1-2270 


10 50-0 


14 13 


n- 7 


1-8426 


1-8633 


8 47-5 


10 20 


n = 9 


2-6815 


2-6951 


7 18-6 


8 7 


?& — X X 


3-7193 


3-7287 


6 12-9 


6 41 



n 


■ 


1 


n 


z=z 


3 


n 


= 


5 


n 


= 


fy 


n 


-~ 


9 


n 


1 


11 



n 


' 


X 


n 


* 


3 


n 


— - 


5 


n 


= 


7 


n 


zzzz 


9 


n 


-, „„„ 


11 



Depth, 14,520 feet (%/4co 3 a 2 = -£$) ; p/cr = '18093. 



•2819L 
1-0201 
1-9132 
3-1919 

4-8672 
6-9414 



•22204 
1-0160 
19300 
3-2025 
4-8741 
6-9461 



25 


28-0 


11 


54-3 


8 


38-3 


6 


423 


5 


26-1 


4 


332 



Depth, 29,040 



leet 



•36381 
1-57362 
3-3391 
5*8906 

9-2387 
13-3856 



(hg/4<ji)*a? 



•32658 
1-57822 
3-34816 
5-8959 
9-2421 
13-3880 



~ To); p/<r = '1*4093. 



20 
9 
6 
4 
3 
3 



599 
33-1 
33-5 

565 
56-8 
16-8 



Depth, 58,080 feet (hg/ico^ = |) ; p/a = -18093. 



n = 1 


-52763 


•50650 


16 51-7 


72- = 3 


2-6806 


2-6853 


7 19-4 


n = 5 


6-1911 


6-1957 


4 49-3 


» = 7 


11-2879 


11-2906 


3 34-3 


w = 9 


17-9816 


17-9833 


2 49-8 


w. = 11 


26-2742 


26-2753 


2 20-5 



41 


55 


16 


8 


10 


3 


7 


18 


5 


44 


4 


44 



29 


39 


11 


25 


7 


7 


5 


10 


4 


4 


3 


21 



20 


58 


8 


4 


5 


1 


3 


39 


2 


52 


2 


22 



§ 9. Numerical Computation of the Height of the Surface- Waves. 

We have next to evaluate the quantities C 3 , C 4 . . . ; when once the periods have 
been determined this will present no difficulty. Suppose w§ are dealing with the 
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type whose frequency approximates to the root of the equation L ;i = ; we have 

seen in the preceding sections how to evaluate K,j + 3 , K // + 4 , . . . and H w _ 3 , H„_ 4J . . , 

Also we have 

C, + 2 /C r = (2r + 3) (2r + 5) K r + 2? 

C,„ 2 /O r = (2r - 3) (2r - 1) H r _ 2 , 



and therefore 



■ * IWNIMH«M> 

vy^ + 2 — 

p — 






(2ti + 3)(2n+ 5)K,,C„ 

(2w + 3)(2rc + 5) (2w + 7) (27i + 9) K„ + 2 K rt + 4 C„ 



• * * * 



o 9 



(2»- 1)(2»-3)H„_ 2 C„, 

(2n — l)(2n — 3) (2w - 5) (2n - 7) H M _ 3 H„_ 4 C„. 



Thus the height of the surface-waves is given by 



£ = C B e at 



" + (2w - 1) (2w - 3) (2n - 5) (2» - 7) H fi _ 2 H,_ 4 P.„_ 4 

+ (2re - 1) (2w - 3) H»_ 3 P B _ 3 + P„ + (2» +3) (2» + 5) K ;t+2 P Ji+2 
+ (2n + 3) (2w + 5) (2w + 7) (2n + 9) K M+S K„ +4 P s+4 + , . . 



where X is the root of the frequency-equation in question, and C„ is an arbitrary 
constant. 

Continuing with the particular numerical example dealt with in § 7, we take 



4^T 3 



= 2'25394, or 



and deduce 



-J.T — 

,<J 



i 



L 
L 6 

La 

1 1 



10 



•281944, 
"065179, 
•018021, 
•000232, 
•008378. 



X 
2w 



= 1-5014 




XtJ 


•013199, 


JLi l J ~— ™ ^^ 


•016171, 


L 16 = 


•01814, 


,Liia ^ "-"" 


•01951, 



Neglecting K 20 , we find 
log K 18 = n 5'535 3 



log K 16 = n 57785, 



log K 13 = n 4*4397, 



and in like manner 



log K H = n 4-0698, 



WK 10 = W4-9812, 



log Ho = 3*2064, 



loff H, 



3-0458 



log H„ =3-0753 
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from these we deduce 

log (C 2 /C 4 ) = 27505, log (C 10 /C 8 ) = n T-5822, 

log (C 4 /C 6 ) = 1-0414, log (C 12 /C 10 ) = « 1-1994, 

log (C 6 /C 8 ) = 1-3653, log (C H /C 12 ) = n 2-9636, 

log (C 16 /C w ) — n 27885, 
log (C 18 /C 16 ) = w 2-647; 



whence, finally 








C 2 /C 8 = -0014, 


CW^s — — "3821, 




C 4 /C 8 = -0255, 


C 13 /O g = + "0605, 




C fi /C 8 = -2319, 


C 14 /C 8 = ~ -0056, 
C 1C /C 8 = + -0003, 
C 18 /C 8 = - -00002. 



Combining with our solution a second, obtained by changing the sign of i where- 
ever it occurs, we obtain a solution in the real form 



£ = C 8 cos (Xt + c) 



•0014P 3 + *0255P 4 + '2319P + P 8 

- -8821P 10 + -0605P I3 - '0056P 14 + -0003P lfl 



where C 8? e are arbitrary constants. 

This determines the type of oscillation for that particular mode which is in question. 
It will be seen that the coefficient of P 8 predominates, and that consequently the 
deformation of the surface will be similar in character to that which takes place when 
there is no rotation, in which case the height of the surface-waves is expressed by a 
single harmonic term. The nodal circles will however be displaced from their 
positions when the rotation is annulled. 



§ 10. Numerical Expressions for the Height of the Surf ace- Waves. 

By the method illustrated in the preceding section I have computed the series 
which indicate the types of oscillation for each of the forty-eight cases for which the 
periods are tabulated in §§ 7, 8 ; these series are given in the following tables. To 
obtain the height of the surface-waves, the series here given must be multiplied by a 
simple harmonic function of the time of arbitrary amplitude and phase, but whose 
period is found from the corresponding entry in the preceding tables. 
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Table III. — Heights of Surface-Waves for Symmetrical Types. 

hg /laPa* = -^ ; (7,260 feet). 



M 


— 


2 


n 


rrr 


4 


w 


zzz 


6 


% 


z^: 


8 


w 


rr: 


10 


7i 


= 


12 









P 3 - 1-2678P 4 + -52G7P 6 -~ -1097P 8 4- -0137P ]0 - '0012P 13 + -0001P 14 - . . . 

2373P 3 + P 4 - -8753P 6 + -2595P 8 - -0403P 10 + *0039P 12 - -0003P 14 + . . . 

0269P 3 + *2714P 4 + P - -5453P 8 + -1139P 10 - -0132?^ + -0010P 14 - -0001P ]6 + . . . 

0014Po + -0255P 4 + -2319P f) + P 8 - -3821P 10 +.-0605P 13 - -0056P 14 + '0003P 16 
. . + •0012P 4 +*0196P 6 + -1977P 8 + Pio- , 2934P 13 +-0373P u --0028P 16 +'0001P 18 
. . +-0009P 6 + -0150P 8 + 4715P 10 +P 12 --2380P 14 + -0252P 16 --0016P 18 + *0001P 20 - . . . 



hg/AaPa? 



YqI (14,520 feet). 



n 
n 
n 
n 
n 
n 



2 
4 
6 
8 

10 
12 



P 9 - -7484P, 4- -1707P* - -0188P R + '0012P in - -OOOIP™ + . . . 



10 



13 



1286P a + P 4 - '4070P 6 + -0594P 8 - -0046P 10 + °0002P 13 — . . . 
0069P 3 + -1311P 4 4- P 6 - -2556P 8 + -0262P 10 - -0015P 13 + -0001-P 14 - . . . 
0002P 3 + -0062P 4 + '1127P 6 4- P 8 ~~ '1840P 10 + «0144P 12 - -0007P 14 + . . . 
. . + -0002P, + -0048P, -|- '0969Po + P 10 - -1432P™ + -0090P U - <0003P lfi + . . . 



. . + -0001P 6 + -0037P 8 + -0845P 10 + P 13 ^ 'lirOP 14 + '0062P 16 - 0002P 18 + . . - 



hg/i 



9, ° 



JL 9 

TO ? 



(29,040 feet). 



n 

n 

n 
n 
n 
n 



2 
4 
6 
8 
10 
12 



P 3 - '4029P 4 + -0477P 6 - -0027P 8 4~ -0001P 10 - . . . 
0677P 3 + P 4 - -1989P 6 + -0144P 8 - -0006P 10 -f . . . 
0017P 3 + *0651P 4 H- P 6 - -1259P 8 4- -0064P 10 - '0002P 13 + . . . 



. . + -0015P. 4~ -0560P- + P« - -0912P 10 + -Q036P 



8 



10 



12 



. . 4- -0012P 6 + -0482P 8 -f P 10 - 
. . 4- -0009P 8 4- -0421P 10 + P 13 



•0712P 13 4~ -0022P 14 - 
- -0583P 14 + '0015P 1G 



•0001P la . + . . . 



hg/AaW = J ■; (58 ? 080 feet), 



n 
n 
n 
n 
n 
n 



2 

4 
6 
8 

10 
12 



P 3 - -2076P 4 4- '0125P G - -0004P 8 + . . . 

0347P 3 4- P 4 ~ '0989P 6 4- '0036P 8 - '0001P 10 4- . . . 

0004P 3 4- •0326P 4 4- P 6 - -0627P 8 4- *0016P 10 - . . . 
. . 4- -0004P 4 4- '0280P 6 4- P 8 ™ -0455P 10 + 0009P 12 - 
. . 4- -00,3P 6 4- -0241P 8 4- P X o - '0355P ]3 4- -0006P 14 - 
. . 4- '0002P 8 4- -0210P 10 4- P 12 - '0291P 14 + -0004P 16 
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Table IV. — Heights of Surface- Waves for Unsymmetrical Types. 

hg/'LcoW = -& ; (7,260 feet). 



n 


— 


1 


n 


m 


3 


n 


— 


5 


n 


— 


7 


n 


rr: 


9 


n 


— 


11 









Pj - 1-5058P 8 + -7106P 5 - -1673P 7 + -0235P 9 - -C022P n + '0001P 13 

1221P X + P 3 - 1'0907P 5 + '3875P 7 - '0701P 9 -f- -0077P n - 

0161P! + '2772P 3 -f P 5 - -6837P 7 + '1688P 9 ~ -0226P n 4- 

OOlOPi + '0280P 8 + -2521P 5 + P 7 - '4498P 9 + -0811P U - 

. . + -0014P 3 +-0224P 5 +-2137P 7 + P 9 --3320P n +'0468P 13 

. . +-0010P 5 + '0171P 7 +*1837P 9 +P 11 ~-2628P 13 4-'0304P 15 



- -0006P 13 + . . . 
•0019P 13 - -0001P 15 + 
•0083P 13 + -0006P 15 - 
-•0039P ]5 +-0002P 17 - 
-•0021P 17 -f-0001P 19 - 



hg/iata* = -£$ ; (14,520 feet). 



n = 


1 


P x _ 1-0477P 3 4- "2986P 5 - *0396P 7 + -0030P 9 + -0001P U 4- . . . 


n = 


3 


•0778P! + P 3 - -5478P 5 + -0993P 7 - -0091P 9 4- -0005P U - . . . 


n = 


5 


•00481^ + -1374P 3 + P 5 ™ -3156P 7 + '0381P 9 - '0025P U 4~ -0001P 13 


n = 


7 


•OOOlPi + -0068P 3 + -1218P 5 + P 7 - '2141P 9 4- '0190P n - '0010P 13 + . . . 


n = 


9 


. . . + -0002P 3 4- '0054P 5 + -1043P 7 + P 9 ~ '1611P U + '0113P i3 - -00Q5P 15 + . . . 


w = 


11 


. . . + -0001P 5 + *0042P 7 + -0903P 9 + P n - -1288P 13 4- '0074P 15 - -0003P-, 7 + . . . 



hgjlaW = iV ; (29,040 feet). 



n = 


1 


Pi - -6516P S + '1034P 5 - -0073P 7 4~ '0003P 9 - . . . 


n = 


3 


•0471P! + P 3 - -2726P 5 4- «0248P 7 - '0011P 9 + . . . 


n = 


5 


•0013P! + -0689P 3 + P 3 - -1547P 7 + -0093P 9 - -0003P n + . . . 


n = 


7 


. . . 4- -001 7P 3 + *0605P 5 4- P 7 - *1059P 9 4- -0047P U - -0001P 13 + . . . 


n = 


9 


. . . -f -0014P 5 + -0519P 7 + P 9 - -0800P n + -0028P 13 - '0001P 15 -f . . . 


n = 


11 


. . . + -0010P 7 -f -0450P 9 + P n - -0641P 13 + -0018P 15 - . . . 



%/4co% 2 



5 5 



(58,080 feet). 




Pi 



. 3 « 



7P, + '0310P, - -0011P 7 + 



0265Pi + P 3 - 1361 P 5 + -0063P 7 - -0001P 9 + . . . 
0003P! + *0346P 3 + P 5 - *0770P 7 + -Q023P 9 - . . . 
. . + -0004P 3 + -0302P 5 + P 7 - *0528P 9 4- -0012P U - 
. . 4- '0003P 5 + «0259P 7 + P 9 - -0399P U + -0007P 13 - 



tat 



. . 4- '0003P 7 + -0225P 9 4- P n - -0320P 13 4- '0005P 
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§11. Forced Tides. 

Leaving now the problem of the free oscillations, let us return to the equations of 
§ 5, when we retain the y$. It is obvious in the first place that a disturbing force 
whose potential at the surface is expressible by surface-harmonics of even order alone, 
or of odd order alone, will give rise to a forced oscillation of like character. Further 
we may consider separately the effects of the different terms in the disturbing 
potential and superpose the results. Suppose for example that the surface- value of 
the disturbing potential is expressible by the single harmonic term 

yj 3 , (/x) e M 
where we will suppose n even. 

The equations (23) which determine the type may be written 



- C^t + CJ7 . 9 = 

cys . r ~~ ql* + cy 1 1.13 = 

0„_ 2 /(2n - 3) (2» - 1) - C„L, + C„ +8 /(2n + 3) (2n + 5) = yJi/AM 
C/ftn + 1) (2n + 3) - C„ +3 L W+S + C w+4 /(2n + 7) (2w + 9) = 



with the condition that C„ = ; whence we obtain 

Cr-t/Cr = (2r - 3) (2r - I) H r _ a (r < w + 1) 

C, +a /0r = (2r + 3) (2r + 5) K,. +2 (r > w - 1) 
and therefore 

C» {H*_ 3 + K fH % — L„J = yji/ioj^, 
or 



C 



«- 



4w% 3 (H^g + Iv w+3 - L n ) 



Thus the height of the tide is given by 



y* h 

4© 2 a 2 (H„„ 3 + K n+3 — L n ) 



; . . + ( 2n - 1) (2w - 3) (2n - 5) (2w - 7) H^H^P^ 
+ (2n-l)(2n-3)H^ 2 P^ a +P l( +(2n + 3)(2w + 5)K, +s P jr+8 
+ (2n + 3) (2n + 5) (2n + 7)(2n + 9) K, +3 K w+ ,P #+4 + . . . 



The expressions H, K, L all depend on X the frequency of the disturbing force. It 
is obvious from the above that the tides become very large when X approaches a root 
of the equation 

&4}i "* OM * . J7|.^ w ****** J\.^4„o »*•«■ 
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and this equation, as we have already seen, is the equation which determines the 

periods of free oscillation. 

It is usual in Tidal Theory to express the height of the forced tides in terms 

of the height of the corresponding " equilibrium-tides." If we denote the height of 

the equilibrium-tide arising from the disturbing potential in question by CU?» (ft) e /A *, 

we see, on omitting all the terms on the left of (23) which depend on inertia, 

and replacing C u by (§,„, that 

hg t fi u j^ar' = yJi/4o)Xt 2 ; 
and therefore 

If then £ denote the height of the equilibrium- tide, we have 



L 

y 

to 



Jig n /4:Gra? 



\ti-u~2 ' -^" + 2 ' J »/ !?■ 



+ (2n~l) (2w-3) (2»-5) (2n-7) H^H^.; 

+ (2n-l) (2«-3) H W _ S P„_ 2 +P„ 

+ (2w+3)(2«+5)K w + 2 P M + 3 

+ (2»+3)(2n+5)(2n+7)(2w+9) K„ + s K. +4 P„ +4 

+ ••• 



The most important practical application of the above theory is the case where the 
disturbing potential involves only a single harmonic term of the second order, and 
the period of the disturbance is long compared with the period of rotation. Thus 
taking X^/ico 2 = "00133, /tc//4w 3 a 2 = 1/40, which corresponds to the case of lunar- 
fortnightly tides in an ocean of depth 7260 feet, we find 



4> 


- '09349, 


L 8 = - -03105, 


L t =- 


- -04745, 


L I0 = - -02886, 


L 6 = - 


- -03561, 


L 13 = — -02764. 



whence, neglecting K 14 , we obtain in succession 



log K 13 = n4'12, log K 10 = ni'432, log K 8 = n4 , 8151 
log K = w3'3055, log K 4 = wS'9992. 



Thus 



L<2 — K 4 



log (C 4 /C 3 ) = nl -7986, 

log (cyc 8 ) 



= - '08351, 
log (C 6 /C 4 ) = nI-4608, 
: nl-033, log (O 12 /C 10 ) 



log (C 8 /C 6 ) 
= »2-88, 



nl-2216, 



C»/6, 



40 g 
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and therefore 

. JL _ -2669P 3 — -1G78P 4 + *Q485P (j - *008lP 8 + '0009P 10 - *GQ0lP I3 + . . . 
In the same manner, when hg/iora 2 = -/q, 



when hg/4cx i a" — yg 

r /p 3 



= -4079P, - -1671P, + -0285P G - '0027P, + '0002P 



•5697P, - M388P.J, + -013lP - -p006P 8 + . . . ; 



and when hg/iofla 

K 



2 1 

— 5? 



i 7208P 2 - -0973P 4 + -0048P 6 - *0001P 8 + . . . 

The lunar-fortnightly declinational tides have been evaluated by Professor Dak win # 
for depths which correspond with the first and third cases given above, the results 
being expressed in series proceeding according to ascending powers of the variable /x. 
If we replace the various powers of jjl by their values in terms of the zonal 
harmonics,t we may deduce the following series from those given by Professor 
Darwin ; when hg/icohfi = -^ we find 



10 



« • » 



jTTS- = "2889P8 - -1755P.J, + '0490P fl - -0079P S + 'OG09P 
while, when %/4ara 3 = y^, 

f == *5969P B - 1385P 4 + -0126P« - *Q0Q6P S + . . . 

The difference between these expansions and those we have given above, is to be 
explained by the fact that we have included in our analysis the effects due to the 
attraction of the water on itself. I have re-computed the lunar-fortnightly tides, 
starting with the assumption that p/cr = 0, and obtained practically identical results 
by the two methods. 

We see then that the effect of the gravitational attraction of the water is to diminish 
the tides, as compared with the equilibrium tides, in the first case by about 8 per 
cent., and in the second by about 5 per cent. 

* ' Encyc. Brit.,' Art. "Tides." § 18. 

f Ppkrpps : i Spherical Harmonics,' p. 27. 
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§ 12. Lunar- Fortnightly Tides in an Ocean of Variable Depth. 

A similar method of treatment may be employed when the depth follows the less 
restricted law, of the form k + I (1 — i$)> made use of in § 5. The numerical com- 

nutation is greatly facilitated in this case when I takes the form — -— where r is 

1 & J r (r + l)g r 

a small integer, since, as we have seen, the series which express the tide-heights will 

then rapidly terminate. 

For example, taking r = 4 so that 



Iff 


1 

~ 4 . 5gjc/ ~ 




1 




i 


4.5 • 


OCT 





which makes the value of / for the earth about 15,454 feet, the values of C 2 , C 4 are 
given by the equations 



* s l4»% a "*" 7 4o)% 2 J 73 UwV- "•" 7 4ft> V 
£ P __ T P ~ — ^ ^ 2 ' ° ^ ft 



K 



whilst all the remaining C s vanish. The notation employed is that introduced at 
the end of § 5. 

ha 1 

Taking j— 3 = — , p/cr = '18093, we find in the case of the lunar-fortnightly tides 

L» = — '27G60, L 4( = - '23787, 
log &= 2-3105, 



whence we obtain at once 



kj , 4 k \ 9% ®s 



C 2 = - /f- + ~-^- ^f 2 ^ '7426(L, 
^ \4ft) 3 a 3 7 4ft) 2 a 3 / # L 3 2 ' 

b * = 577 w lI + It 3 = " ' 0980( ^ 

iXiCJ 1 I ft 

Thus, if j-y-- = — -f 7-^ — sin 3 0, where 6 denotes the co-latitude, which for a 
system of the dimensions of the earth makes the law of depth 

(58080 + 15454 sin 3 6) feet, 
we find for the tide-height the expression 

£= (S a [7426P a - -0980PJ. 

MDCCCXCVII. — A f 2 I 
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Similarly, when 



when 



a 



nd when 



7 4o)% 3 / 1 1 ci 

h = - -r + — - — 

// \10 r 4 . 5 ^ 



sin 3 0) = (29040 + 15454 sin 2 0\ 



1= 6 3 [-620lP a - '1446P 



4 



/i_ g Uo + 4.5^ &m " 



g = 6 8 [-5018P 3 - -1897P, 



(14520 + 15454 sin 3 0). 



h = 4 — f r« + -— - 7 siu 3 ) = (7260 + 15454 sin 3 0). 
g \40 4.5# 4 



Again, if we take 



J = g a [-4095P 3 - '2248P t 



A 2 2 



the series for £ will terminate with a term involving P 6 , and the values of C 3> C^, C 6 
must be computed from the equations 



LiA + %£ 



m 



4gA& 3 ""*" 7 4o)% 2 / </ — \4o>% 3 "^ 7 IgA* 8 / 2? 



^2^2 ~~ -1^4 



b-i^i- ^0(5 



2 . 3 ly% 

577 4«% 3 



2 



3 4/3 



5,7 4ct>% 



o o ^2' 



= • 



from which we find, when j^j^ = -y , 



4&) 3 a ! 



= '7316(5, 



2? 



G 



4 



•0978(5 



25 



C fl = '0024(5 



'6 



2< 



Thus for the law of depth 



4o)% 3 f 1 , 1 fj . 



5 



// L o 6.7 # 6 
g = (5 3 [7316P 2 - -0978P 4 + *0024P 6 



sm 3 0!> = 58080 + 7218 sin 3 0, 



in like manner, when 



h 



4a>% 



9.„. a r 



9 



1 



1 



- J- -■ — sin 3 
10 + 6.7 /7. 



29040 + 7218 sin 3 0, 



£ = (S 8 [-5954P S - -i426P 4 + -0064P fi ] ; 
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when 



and when 



h 



^coht 2 r i 



9 



H~ n 



1 9 • 



20 ' 6.7 g fl 



sin 3 6 \ = 14520 + 7218 sin 3 0, 



I = g 2 [-4576P s - -1821P, + -0137P fi ] ; 



h 



4:C0 2 d 



m r 



i *- // •0/1 

1 7a + ^ ' hr — sir 
[40 6.7 ( h 



= 7260 + 7218 sin 2 0, 



1= 6 3 [-3457P a - -2075P 4 + -0234P 6 ]. 



For other values of I we must employ a method similar to that of the last section. 
The general formulae for the computation of the forced tides due to a disturbing 
potential of the second order are 



where 



L A + 77 2 c 4 = i7^~ + 



4- % 



&/, 



4o)% 3 ' 7 4«V 
f 4 C 4 - L 6 C fl + >? 6 C 8 = 



+ 



4 lg u 



*\ 



9 -.9 



4w 2 a 2 ' 7 4w% 



(So 



2.3 



kh 



5 . 7 4co 2 « 



o6 2 



> 



J 



• • 



<J» 



1 — n (n + 1) lg n /4:G) 2 a? 
(2n + l^(2n + 3) 



^-2 



— 1 - ^ (^ + jO IgnfiaW 



L w — 



(2rc- 1)(2t& + 1) 

i_ i 2{ l — 7t( ^+ l) lg n /4:(D 2 a ?} kg n 

n(n-t 1) + (2rc ~nL^(2?7+~;3J 4^>V 



(31), 



Let us introduce the notation 



JUL 



% ~ 



K 



%'iVn 


§ 


fc 


»-2^»-3 


L« 


— 




Ju»_2 


ft 


S^»" 


-3 


£nVn 






W- 



L 



■H 



L M . n — 



'» + 3 



Then it may be shewn, as in the last sections, that for values of n greater than 2, 



^tf + 2/ ^« — ■"*« + 2/v^ • • 



* ft • * V 



(32), 



and therefore the first two of equations (31) may be written 



f hg 2 , 4 Z#, 



LA + %C 4 = {^ + y 4w i} «a, 



2 l @o. 



^oV,^ 



rtW I. fe 
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On solving these equations we obtain 



a 






kh 1 % 



4am 



9 9 



i. 



7 4ft)V- L, -K 



fr-|-*im 



~d< 



c. 



4 



1 



kh , 4 



^2 1 K 4 S 2 



^.»J Wo -^-2 


6* 




5.7 4o)% 3 £ 3 L 4 - 


-Ho- 


"K 8 


2.o % 


% 





. (33), 



4»% 3 



7 4o> 8 a» J % L, - K 4 5.7 4» 9 a* L 4 ~ H 2 - K 



6->* 



after which we can deduce C 6 , C 8 . . .in succession by means of (32). 

For example, taking 

kg 1 Ig 1 

~4«**a? = 20 ' 4eo 3 a 3 ~ ~ ' 



30 



which makes the law of depth for the earth 



(14520 + 9680 sin® (?) feet. 



we find for the lunar-fortnightly tide 



L 3 = - -13276, L 4 

L 8 = - -07156, L 10 



•08722, L 6 = — -07583 ; 
•06955, L 12 = — '0685, 



and 



log &%■ = '1-1437, log ^ = n 6-9549, log &% = 6-9586, 
log i 8 ij 8 = 5-421 7, log £ l(Wl0 = 5-578, 

whence, if we suppose K 14 = 0, we obtain in succession, from the formula 



K = 



-11 



?»--2 Vn~~J 



»+2 



log R 13 
log K c 



w4'75, log K 10 = n 4-583, log K 8 = n 4-1062, 
4-0758, log K 4 = w 3-2025. 



Also 



log H 8 = n 3*0206, log £,= 2-3707, log ty, = 37730 ; 



whence from (33) we find 



C< 



'4719(5,, C. 



2> ^4 



•• 18526 



2' 



gain 



log(C fl /C,) =log(K 6 /^) =« 2-6976, 
log(C 8 /C 8 ) =log(K 8 /, 8 ) = 2-3910, 
log (C 10 /G 8 ) = log (K 10 / % ) = 2-774, 
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whence finally 

C =:'0092, C 8 = *00023 ? C 10 = '00001 ; 

and the expression for the tide-height becomes 

£ = S a [-4719P S - '1852P^ + *0092P 6 + -00023.P 8 + *0000lP 10 + . . .]. 

As a further illustration, I have computed the series for the tide-height for the 

4ft) 3 ^ 3 / 1 1 . \ . 

case where h = ( ^ — ^ sin 2 Ou that is, where the depth is 14,520 feet at the 

poles and shallows to 4840 feet at the equator. The value of £ in this case is as 
follows : — • 

{ = 6 3 [\3082P 2 - '1 106P 4 + *0467P 6 - '0158P 8 + '0048P 10 

- -0014P 12 + '0004P 14( - -0001P 16 + ...]• 

When I is positive, that is, when the depth at the equator exceeds that at the 
poles, the series appear to converge more rapidly than when the depth is uniform, 
but the opposite is the case when the water is deeper at the poles than at the 
equator. 

§ 13. Forced Oscillations of Infinitely Long Period* 

If we suppose X so small that we may neglect \ 2 /4g) 2 , we find, on putting X a /4(w a == 0, 
for the height of the forced tides the following four series in place of those given 
in S 11 : — 

v %^ = '2661P 2 - -1671P 4 + '0482P fl - *0080P 8 + -0009P 10 - -000lP 12 + . . . 

-|r = '4070P 3 - -1666P 4 + -0284P fi - '0026P 8 + -0002P 10 - . . . 

L. _ . 56 89P 3 - -1385P 4 + -0130P - '0006P 8 + . . . 

So/* 2 

A = .' 7 201P 3 - •0973P 4 + '0048P 6 - '0001P 8 + . « . 

So/*9 

The lunar-fortnightly tides therefore differ only very slightly from tides whose 
period is infinitely long, The difference between these latter and the solar semi- 

* Several of the conclusions of the present section have been previously arrived at by Professor 
Lamb (* Hydrodynamics/ chapter viii.) ; ' but, on account of the important light which they throw on the 
later sections, I have thought it desirable to treat the questions in some detail, even at the risk o 
repeating what is already well known. 
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annual tides will be quite inappreciable, and we may take the above series as giving 
a good representation of the solar long-period tides, unless the effects of friction 
become important for such tides. 

The fact that when the period of the disturbing force is increased without limit 
the free surface does not tend to approach its equilibrium form appears at first 
sight to be at variance with the general laws of oscillating systems. The explanation 
of this apparent anomaly may perhaps be made clear by considering a simple form 
of "gyrostatie" system which possesses only two degrees of freedom. In the absence 
of frictional forces, the general equations of motion of such a system may, by a proper 
choice of coordinates, be expressed in the form 



x ■— o)y + n % x = 

y + ^ x 4™ w % ^ Y. 4 ' 5 " 

Here x, y denote the generalized coordinates of the system, Of the terms on the 

• • • • ® « 

left, the terms x 9 y are due to inertia, the terms cot/, cox are described by Thomson 
and Tait as "motional" forces, and the terms n*x, m 2 y as "positional" forces; X, Y 
are the generalized components of the external disturbing force. 

If now x, y, X, Y be supposed proportional to & x \ we find from the above equations 

— k^x — <oiky + n^x = X, 

— k 2 y + ankx + m z y = Y, 

whence we may obtain x, y in terms of X, Y» When the period of vibration is 
indefinitely prolonged, k will approach zero as a limit, and the limiting form of the 
solution will in general be 

x = 'X/ri 2 , y = Y/m 2 . 

This implies that the displacements will in general tend to acquire their equilibrium- 
values as the period of the disturbing force is lengthened. There will, however, be 
an exception to this law if one or both of the positional forces n 2 x, m 2 y vanish. 

Let us first examine the nature of the free oscillations in such cases ; omitting 
X, Y, and supposing that n = while m remains finite, we have for the determina- 
tion of the free motions 

— k?x — mky = 

— k 2 y + (oikx + whj = 



> 



or, if we denote by u, v the generalized velocity-components so that u = x = ikx. 

v = y = iky, 

* Thomson and Tait, 'Natural Philosophy,' vol, 1, p. 396 (1886 edition), 
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i\ (u — coy) = 0, 
— K 2 y + cou + wify = 0. 

'Til 'II 

These equations will be satisfied if A. = 0, u- = = const. lb follows that the 

system is capable of a small free steady motion relative to the rotating axes, defined 

by 

y =. const, it = —■ mhjjco. 

If both n*x, m l y are zero, the equations for the free motions become 

— X 2 x — (oiky = 0, 

— X 2 y + coiXoj = ; 

both of which are satisfied by supposing that x, y are small arbitrary constants, and 
therefore X = 0. 

In the latter case the equilibrium-state defined by x = 0, y = is not the only 
condition of relative equilibrium, but any other configuration of the system in the 
neighbourhood of this one will also form a configuration of relative equilibrium. 

Let us now consider the nature of the limiting forms of the forced oscillations 
when the period is indefinitely prolonged. In the former case we must suppose that 
the disturbing forces are such that they do no work when the coordinate x is varied, 
so that X = 0, as otherwise the stability of the system will be destroyed ; the 
equations of motion for the forced oscillations then become 



whence 



and in the limit 



iXu — 


onXy = 0, 


-X 3 */ + 


ou + my = Y ; 


u = coy - 


Y(jO 


Yco 


Y 


m* -f- 


u — — « 



The velocity-component u will therefore always remain of the same order as the 
disturbing force Y, while the amplitude of vibration of the coordinate x will tend to 
increase without limit. 

In the latter case the equations of disturbed motion may be written 

i\u — cov = X 
iXv + ton = Y 
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whence in the limiting case 



U = Y/tO, V = — Xy 



CO. 



Hence both the velocity-components tend to finite limits while the amplitudes of 
vibration of both coordinates increase without limit. 

The essential characteristics of both cases are (i) that one or more of the generalized 
coordinates does not appear explicitly in the equations of motion but only the 
corresponding velocity-component, and hence (ii) that X = is one root of the 
frequency-equation for the free modes of vibration, from which it follows that 
(a) either free steady motions relatively to the rotating system are possible, or 
(/3) that the configuration of relative equilibrium defined by x = 0, y = is not 
isolated. The two conditions (a), (/J) may both be expressed by stating that the 
steady motion defined by x = 0, y = is not the only form of steady motion of which 
the system is capable. 

The two cases are both illustrated by our problem. For if we suppose the waters 
of the ocean displaced horizontally in such a manner that the form of the surface is 
unaltered, we shall evidently obtain a new configuration of relative equilibrium, while, 
as we shall see in the next section, if we suppose that the fluid is in relative motion 
in such a manner that the fluid particles are moving along parallels of latitude, it is 
possible by a proper adjustment of the free surface to ensure that such a motion should 
be permanent. 

The coordinates which depend on the horizontal displacements alone are analogous 
to the coordinate x in the former of the illustrations we have given above, and to the 
coordinates x, y in the latter. They do not appear explicitly in the equations of 
motion, but only through the corresponding velocity-components. We conclude, then, 
that the horizontal velocities will be ov the order of the disturbing forces, whereas the 
horizontal excursions of the fluid particles will tend to increase without limit as the 
period is prolonged. 

By way of explaining how these circumstances may arise physically, let us suppose 
for the moment that X is actually zero, and consequently that the disturbing force is 
constant. In the case of a system oscillating about a position of equilibrium, the 
introduction of a constant disturbing force will have the effect of slightly changing 
the configuration about which oscillations corresponding with the free modes of 
vibration take place. Suppose now that a disturbing force, such as that which gives 
rise to the long-period tides, tending to increase the surface-ellipticity of the ocean, 
is suddenly applied to our rotating system when in a configuration of relative equi- 
librium. It will immediately set up oscillations, the initial motion being such that 
each particle will tend towards the position in which it would be in equilibrium under 
the new disturbing influence. The new position of equilibrium is such that in it there 
will be more water in equatorial regions, and less water in polar regions, than in the 
old. Thus the initial motion involves a flow of water directed from the poles towards 
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the equator. The water however coming from higher latitudes into lower will reach 
these lower latitudes with an amount of rotation less than that which is appropriate 
for these latitudes if the whole were in a state of steady motion as a rigid body. 
There are no forces acting' which tend to modify the angular momentum about the 
polar axis of an elementary ring of water, which coincides with a parallel of latitude, 
and consequently currents will be started, in virtue of which each particle of fluid 
will move along a parallel of latitude from east to west. The effect of the disturbing 
force is therefore to modify the state of steady motion about which the free oscilla- 
tions take place from a uniform rotation of the whole system as a rigid body to a state 
in which there exist horizontal westerly currents. If, as is usual in dealing with 
forced oscillations, we suppose the free oscillations to be annulled, we see that the 
" forced oscillation " arising from such a constant disturbance as we have been con- 
sidering will be of the nature of a steady motion relatively to the rotating earth, 
consisting of a westerly flow of the whole ocean, the velocity however varying with 
the latitude. 

In the case of a periodic disturbance of very long period, the motion set up at any 
instant will be of like character, provided that the viscosity of the fluid is not 
sufficient to sensibly affect the currents in question in the course of a single period. 
An equilibrium-theory will only be applicable when the rate of dissipation of such 
motions is so rapid that they practically disappear in a time which is short compared 
with the period of the disturbing force. Now in an ocean whose depth is equal 
to the mean depth of the actual ocean, it seems highly improbable that such 
currents would be appreciably affected by viscosity in the course of a few months. 
Hence it appears that the present theory in which the effects of viscosity are totally 
disregarded will almost certainly give a far better representation of the lunar long- 
period tides than the equilibrium theory, and most probably also of the long-period 
solar tides. 

§ 14. Free Steady Motions. 

In the last section we have called attention to the fact that free oscillations of 
infinite period are possible, or that the system with which we are concerned is 
capable of free steady motions. We proceed in the present section to examine the 
nature of these steady motions. 

Referring back to §§ 2-4, we see that the general equations of motion of the ocean 
when free from external disturbing influence, and at the same time supposed steady, 
so that \ = 0, can be expressed in the form 
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dt ~ 

u = 

V = 



a 



i 

d/ju 



{/(l-^)AU) + ~( 



hV 



d<j> W (!-/*») J J 



-\ 



2ft)«/^ y (1 — /a 3 ) 9<£ 

2o>^ 3yL6 ? 



1> 



V 



rf 



J> 



(34). 



where, as there is no longer any ambiguity f we have omitted the bars from the 
symbols U, V, \jj, v'. 

Since 3£/3£ = when the motion is steady ? we see that the first equation is 
identically satisfied if we suppose h 9 \jj both independent of <j). In this case we may 
take 



where the constants C» are arbitrary, and deduce 



Q 9 9 » * 



* » * • \^*^/5 



whence 



$= —t gnOn'Pu 0*) ; 






. (36). 



The last equation gives the velocity which must be imposed on the particles of 
water in latitude sin"" 1 /*, in order that the free surface may be maintained in the form 
defined by (35) without any external force. We see that it is theoretically possible 
to maintain an arbitrary surface-form by correctly distributing the longitudinal 
velocities of the fluid particles. If however the series (35) involves harmonics of 
odd order the value of V given by (36) becomes infinite at the equator, and to 
prevent a flow of liquid across the equator it would be necessary to impose an 
infinite velocity on the particles of water there. Hence, if the water extend either 
wholly or partially over both hemispheres, the distribution of velocity and the form 
of free surface must be symmetrical with respect to the equator, at least so far as 
concerns that part of the ocean which communicates across the equator. 

Conversely, any arbitrary initial distribution of longitudinal velocity symmetrical 
with respect to the equator may be rendered permanent by an appropriate adjust- 
ment of the free surface. These results hold good whatsoever be the law of depth, 
provided it be a function of the latitude alone. 

If \]j be not supposed independent of <f>, we see from the second and third of 
equations (34) that the velocity of flow across any element ds inclined at an angle x 
to the meridian is 



ANALYSIS TO THE DYNAMICAL THEORY 0¥ THE TIDES. 



251 



U sin x — • V cos x 

1 



- «/ 



2coa/jL 

1 3^ 
2oy/M ds * 



sm X 



1 djr 



+ cos x \/(l — /^ 2 ) 



3i/r 
3//, 



Hence, if 3t///3s = 0, there will be no flow across the element ds. It follows that 
the function xjj may be regarded as a stream-function, the paths of the particles of 
water always coinciding with the lines 

i/f = const. 

But from the equation of continuity we have, on putting d^/dt = 0, and replacing 
U, V by their values in terms of \jj, 



= 



¥ 



lb 3\Jr 1 3 [ h 3-*|r 

/ju 3<£ J d(J> \ fM d/ju 



or 



A / A\ 8 t _ A / A\ a t _ o 



the general solution of which is 



\Jj =f (Jl/[l) . . . (^^): 



where f(h/ii) denotes an arbitrary function of h/p. 

It follows that the stream-lines \ft = const coincide in direction with the lines 

ft J Lit — — COUot/ »»»•«••»•• 



(38). 



Thus, if the depth be a function of the latitude alone, the stream -lines must necessarily 
coincide with the parallels of latitude, and the only forms of steady motion possible 
are those in which the water has no latitudinal velocity. In the more general case, 
the stream-lines of the possible steady motions are given by the equation (38), and 
from this equation they might at once be traced out on a chart if we had a sufficient 
knowledge of the depth of the ocean in different parts. In particular, whatever be 
the law of depth, the equator will be one of the free stream-lines corresponding to an 
infinite value of hjfx, while the shores will also be stream-lines corresponding to zero 
values of this expression. An infinite number of stream-lines will converge towards 
those points where the coast-line intersects the equator, and it is only by passing 
through one of these points that a particle of water could pass from the northern to 
the southern hemisphere, or vice versd. As however the velocities at these particular 
points tend to become infinite 5 the equations which we have used which involve the 
neglect of the squares of the velocities will not be applicable to the region immediately 

2 K 2 
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surrounding the points in question. It seems most probable that the inclusion of the 
terms involving the squares of velocities would have the effect of diverting the stream- 
lines, so as to cause them to follow the coast-lines even in the immediate neighbour- 
hood of the equator. 

An attempt to trace out the lines h/fi = const for the North Atlantic Ocean from 
data obtainable from the Admiralty charts quickly showed that the theory here put 
forward is inadequate in itself to explain even the more salient features of the 
circulation in the region in question. Observation however indicates the direction 
in which we must look for the defects of this theory, The excessively low tempera- 
ture of the water beneath the surface* in equatorial regions can only be explained 
by supposing that this water has travelled thither from higher latitudes, whereas we 
know that the currents at the surface, for the most part, set from the equator towards 
the poles. We conclude that the under-currents of the actual ocean differ materially 
from the surface-currents, and in this respect the actual circulation differs from the 
types of circulation with which we have hitherto been concerned and which are the 
only possible types of circulation in our ideal ocean in which the density is uniform. 
It has been urged by some authorities that the variations in the density of the water 
arising from differences of temperature, salinity, &c, are the sole causes which 
maintain ocean currents, but in that we have seen that currents could exist even 
without such variations, it seems to me to be highly improbable that such is the 
case, though there can be no doubt, in the light of our present analysis, that these 
variations are largely effective in determining the course which the currents pursue. 

If we suppose that the ocean consists of a number of horizontal layers of different 
densities, but that the density throughout each stratum is uniform, then for each of 
the strata a function \jj will exist defined by 

\js = V' + \o>\x % + y 2 ) — — + const, 

P 

and the horizontal velocities for any stratum will be connected with the corresponding 
function by the equations (34). The equation of continuity for any stratum may be 
formed as in § 3, provided we replace h by the depth of the stratum in question 
instead of the depth of the whole ocean. The stream-lines for any stratum will 
therefore still be given by the equation (38), with this modification in the meaning 
of the symbol h. It follows that the equator will still be one of the free stream- 
lines, but the motion elsewhere may be totally different from what it would be if the 
density were the same throughout. It seems probable then that the result we have 
obtained with reference to the tendency of the currents to set alonej the equator will 
still hold good even when the density is variable, and this conclusion is borne out by 

m The principal facts at present known in relation to the distribution of temperature in the ocean 
will be found in the ' Eeport of the .Challenger Scientific Results/ s Chemistry and Physics/ vol* 1. 
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observation, it being noticeable at a glance at the Admiralty current charts that 
there is no tendency to cross the equator except in the immediate neighbourhood of 
the coasts. 

The rigorous treatment of the problem of ocean currents, as affected by variations 
in the density of the water, appears to be hopelessly beyond the powers of mathe- 
matical analysis, and I will therefore leave the subject with the brief indications 
already given in this section, and will conclude the paper with an example illus- 
trating another means by which possibly ocean currents are in part maintained, and 
which is instructive in showing the very important part played by the rotation of 
the earth in rendering effective a cause which otherwise could give rise to no 
sensible currents. 

§ 15. On Currents due to Evaporation and Precipitation. 

A cause which has been advocated # in explanation of ocean currents is the 
fact that in equatorial regions the amount of water evaporated into the atmosphere 
largely exceeds that precipitated in the form of rain in these regions. The excess of 
water in the atmosphere is carried away to be precipitated in temperate and polar 
regions, thereby giving rise to an excess of precipitation over evaporation in the 
latter regions. It has been urged with some reason that, as the actual amount 
of water in equatorial regions does not diminish nor that in polar regions increase 
from year to year, there must be a continual flow of water from the poles towards 
the equator. The fact that this flow of water is in the opposite direction to 
that observed at the surface, which for the most part sets from the equator towards 
the poles, is explained by attributing the counterflow to undercurrents. If however 
we subject the question to the test of mathematical analysis, we shall find that 
though such a flow towards the equator must necessarily exist, it is so slow as to be 
completely masked by larger currents due to other causes. The flow in question 
will however give rise indirectly, in consequence of the rotation, to currents which 
in the absence of dissipative forces would tend to increase without limit. The 
explanation of this fact will be obvious after the discussions of § 11. 

The effects of evaporation and precipitation may be conveniently represented 
mathematically by an appropriate distribution of sources and sinks over the free 
surface. This will modify the surface-conditions at the free surface but will not 

interfere with the dynamical equations. Instead of equating W to d^/dt we must 
replace it by a certain function of the position on the surface, independent of 
the time, but depending on the rate of evaporation and precipitation at the place, 

# See an Article by Proctor in * St. Paul's Magazine,' Sept,, 1869, reprinted in ' Light Science* 
(1st series), p. 114. 
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which over sufficiently long intervals of time we may regard as uniform. A simple 
law which will serve for purposes of illustration may be chosen as follows :— 



W 



a? % (p). 



The equations with which we have to deal will then be 
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(39). 



where, if we neglect the attraction due to the surface-inequalities, we may take 



* 



gl 



To obtain a particular solution of these equations suppose 



* 



9t=-g(Co + CA TJ=TS + V 1 t,V=Y 9 + V 1 t, 



where f , £ l5 &c, are all independent of £. 

Substituting these expressions in the equations (39), and equating coefficients of t. 

we find 
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If /j- be a function of /a alone, these will be satisfied by 



U 1 = 0, V x 






provided £ x be also a function of /x alone. 

Next, if we equate the terms independent of t in the two members of (39), we 
obtain 
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Suppose the system starts from rest in its position of relative equilibrium, so that 
£ = when t = 0, or £q = 0. Then 

U l + 2cofxY = 
V x - 2a> / xU = 0, 



whence, 



v A = El — o 



u. 



^_Yl = _ 9 v/(i-^) er 1? 



and therefore 



(J u 

aP (a) = — -~~« T" 



/x 3 3/i 



The last equation gives 

(1 ~ u 3 ) A 3& 4ft) 2 a 3 f-r* / x 7 

a I P 2 (/*) ^ 
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no arbitrary constant being added, since both sides vanish when ft = ± 1. Thus, 



9£ 2g) 3 & 2 o 
fy 9 H " 
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Suppose for example that h is constant ; we shall then obtain 



coPaPet 



li = I —j^ /* 4 + const 



s) 3 a 3 « 



= i ~Y- (i/V-P* + f -Pa + 1} + const. 

Choosing the constant so that the mean value of £ x over the surface is zero, we 
obtain finally 

/ — 1 C ^L^ f 8. P I 4 T> > 
bl — 2 gfo 13 5 r i* T 7 * 2J. 

Hence the particular solutions of the differential equations which represent the 
" forced" motion due to the disturbing influence in question are 
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2(*y j a?at 



£ = ^t = ~~7^~j~ {2P 4 +* 5P 2 }, 



U = U + U^ = -- \ y ap y/(l — ^\ 

Y - V + V l t= - ^ p* y(l - ^). 

We see then that the effects of evaporation and precipitation will be to cause a 
steady flow of water, not by means of undercurrents only, but by currents sensibly 
uniform throughout the depth, towards the equator ; in addition to these currents 
the cause in question will give rise to longitudinal currents, not of a steady character, 
but increasing uniformly with the time, and these will be accompanied by an appro- 
priate continuous deformation of the free surface. Were it not for viscosity these 
currents would increase without limit and ultimately endanger the stability of the 
system, but under the action of dissipative forces a steady state must ultimately be 
attained, in which the rate at which the currents are generated exactly balances that 
at which they are destroyed. Thus, suppose the type of motion set up is such that 
if left to itself it would be reduced in the ratio 1 : e in a period r. If U denote 
the velocity of any particle, the law of variation of U under the influence of viscous 

forces is then 

dXJ/dt + U/r = S 

whereas, if there be no viscosity and the system is subjected to such a disturbance as 
we have been dealing with, the velocity varies according to the law 

where f is constant. Equating the rate of increase of the velocity without viscosity 
to the rate of decrease under the influence of dissipative force, we find that the 

ultimate state is defined by 

t = /, or U = j To 

Thus, if the disturbing influence tends to set up one of the possible types of motion 
of which the system is capable under viscosity, the ultimate velocity of any particle 
will be that which it would acquire in a period equal to the modulus of decay of the 
type of motion in question. 

By way of numerical illustration, take a year as the unit of time and an inch as 
the unit of length, and suppose a = 40. This will imply an annual rainfall at the 
poles which exceeds evaporation by 40 inches, and an annual rainfall at the equator 
which is less than evaporation by 20 inches. Farther, suppose h/a = -q-^-q-. Then 

U^-iX 2890 X 40 /I y (1 — ja 3 ). 

U will be numerically greatest when fi % = ^, or in latitude 45°, and the 
greatest value of U is 28,900. The maximum latitudinal velocity will therefore 
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amount to 28,900 inches, or about half-a-mile, per annum. This velocity will be quite 
inappreciable to observation. 

The amount of longitudinal velocity generated in the course of a year in any 
latitude is given by the formula 



acta) 



/x 2 v/(l — /x 2 ). 



This will be greatest when y^ = § , that is, in latitude 55° nearly, and its greatest 
value corresponds to a velocity of about \ of a mile per hour. The maximum 
current velocity due to this cause would therefore amount to about four miles per 
hour if the modulus of decay is as long as 20 years. 

The most crucial test to which we can subject the theory of ocean currents here 
put forward will consist in the evaluation of the moduli of decay for the types ol 
motion concerned. This I have endeavoured to do, but as the work involves 
analytical considerations of a somewhat different character from those which occur 
in the present work, I have deemed it advisable to present the results in a separate 
paper. * These results, so far as they are applicable, seem to point to a modulus of 
decay far in excess of the 20 years here required, but the mathematical difficulties 
have compelled me, in dealing with friction, to leave the rotation entirely out of 
account. It appears that in the simpler system so treated, the types of free current 
motion are far more arbitrary in character than those at which we have arrived by 
including the rotation. This arises from the fact that our rotating system will be 
capable of a large number of free oscillatory motions besides those which we have 
examined in which the periods of oscillation always bear a finite ratio to the period 
of rotation. As the period of rotation is lengthened, the period of each of these 
types of oscillation is prolonged, and the possible forms of steady motion where 
there is no rotation must Include the limiting forms of each of these types of 
oscillatory motion. 

The moduli of decay of the free current motions when there is rotation may 
therefore be very different in value from those obtained in the paper referred to, but 
it does not seem to me that they could be much less in order of magnitude than the 
moduli of decay of the principal types of free oscillation. If this prove to be the case, 
the estimate of 20 years, which we have taken for the modulus of decay, will not be 
so excessive as might at first sight appear. 

Of course, if the water does not cover the whole earth, or if the depth be not 
uniform, the currents due to the rotation will follow the free stream-lines defined by 
equation (38) instead of following the parallels of latitude. The rotation will, no 
doubt, produce its maximum effect when the stream-lines in question coincide with, 
the parallels of latitude, but this circumstance does not alter our main conclusion as 
to the adequacy of evaporation and other such causes to generate currents quit 
comparable with those known to exist in the ocean. 

# "Read before the London Mathematical Society, December 10th, 1896. 
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